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INTRODUCTION 


Over  the  last  decade  there  has  been  only  a few  published  studied  on 
the  feasibility  of  using  seismology  for  cave  detection.  Though  the  lit** 
erature  is  limited,  the  results  have  varied  widely.  Cook  (1964)  postu- 
lated that  S-waves  should  strongly  reflect  off  of  large,  brine  cavrtiva, 
since  Zoeppriti  equations  predict  a total  energy  return  for  shear  waves 
reflected  from  rock-brine  boundaries.  He  succeeded  in  developing  a tech- 
nique for  generating  good  S-waves,  but  the  quality  of  the  reflections 
from  the  cavities  was  poor.  However,  he  did  find  amplitude  attenuation 
for  the  portion  of  the  P-wave  which  had  passed  through  the  cavity  in  com- 
parison to  the  portion  that  had  missed  it.  This  anomaly  was  called  the 
seismic  amplitude  "shadow". 

Watkins,  Godson,  and  Wstson  (196?)  investigated  three  phenomena 
possibly  associated  with  near-surface  cavities;  free  oscillations  of 
cave  walls,  seismic  amplitude  "shadow"  cones,  and  anomalous  delays  in 
arrival  times.  A reflection  traverse  over  a lava  tunnel  in  the  Pisgah 
lava  flow  produced  long-lived  oscillations.  Because  they  vibrated  at  a 
set  frequency,  they  were  referred  to  as  cavi'.y  resonance.  Watkins  et. 
a/,  contended  that  this  resonance  was  in  agreement  with  the  radial  os- 
cillations predicted  by  Biot  (1952)  for  a cylindrical  bore  in  an  elastic 
sulld.  Above  the  Kana-a  lava  flow,  they  discovered  some  clear  amplitude 
anomalies  in  the  "shadow"  cone,  but  cavity  resonance  was  not  well-defined. 
At  the  Nevada  test  site,  the  underground  nuclear  nuclear  explosion  U4B 
had  produced  a rubble -filled  chimney  whose  rcuf  was  only  30  meters  below 
the  surface.  The  volume  of  shattered  and  fractured  rock  extended  to  a 
dep  h of  almost  half  a kilometer.  Again  they  recorded  some  amplitude 
anomalies  and  sotr*  possible  resonance.  In  addition,  they  found  delays 
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In  arrival  titles  fot  the  portion  of  the  wave  that  hint  propagated  through 
the  chimney.  Godson  and  Watkins  (1968)  attempted  to  use  cavity  reso- 
nance to  map  a cave  system  causing  severe  leaks  beneath  Anchor  Reservoir 
in  Wyoming.  They  recorded  some  continuous  vibrations  thought  to  be  res- 
onance, but  they  did  not  report  any  success  In  correlating  these  with 
cavities. 

The  field  work  done  by  Rechtien  for  this  pi  eject  consisted  of  seis- 
mic l:nes  laid  across  several  caves  in  east  central  Missouri.  Mo  peri- 
odic, long-lived  signals  t?t&t  could  be  described  as  resonance  were  re- 
corded. However,  some  possible  reflections  from  the  roof  of  a large 
chamber  within  the  Cathedral  Caverns  complex  were  detected.  There  was 
also  some  Horizontal  particle  motion  that  appeared  co  he  polarized  par- 
allel to  the  sides  of  this  chamber  and  its  single  access  tunnel  (Rechtien 
and  Stewart,  1975).  Additional  studies  in  this  area  discovered  some: 
apparent  time  delays  within  the  seismic  "shadow1*  zone  for  sigualn  re- 
flected from  a deeper  horizon. 

The  problems  that  the  theoretical  portion  of  this  research  pro- 
ject attempted  to  solve  were:  (I)  why  should  cavity  resonance  be  sc  dis- 
tinct at  one  location  and  weak  or  nonexistent  at  another;  (2)  what  are 
the  reasons  for  the  anomalous  amplitude  and  time  delays  for  sigralu  re- 
corded above  some  caves;  and  (3)  how  can  the  seismic  signature  of  i cav- 
ity be  enhanced  and  sore  easily  recognized  on  a seismogram?  Several  the- 
oretical approaches  were  employed  in  an  attempt  to  answer  these  ques- 
tions. In  solving  the  wave  equation  analytically  in  various  coordinate 
systems,  Stewart  (1972)  was  able  to  examine  the  theoretical  aspects  of 
radial  oscillations  of  tha  cavity  and  the  interaction  of  this  resonance 
with  the  uarth'a  surface.  Cavanaugh,  Stewart,  and  Rechtien  (1975) 


const Jerecl  several  types  of  cavity  surface  vr av«t  as  being;  the  possible 
source,  of  resonance  found  above  some  caves.  Fle.iural  waves,  which  would 
be  formed  if  the  cavity's  roof  behaved  as  a vibrating,  thin  plate,  were 
also  investigated  as  a possible  source.  Finally,  wave  sea*,  ter  vaa 
thought  to  be  the  cause  of  the  anomalies  found  in  the  "sh  idow"  zone.  To 
substantiate  this  hypothesis,  both  analytical  and  numerical  solutions  to 
the  scatter  problem  wore  ^.udied. 


VIBRATIONAL  MODES  OF  SPHERICAL  CAVITIES 


The  first  systematic  attempt  to  theoretically  investigate  the  pos- 
sibility of  cavity  resonance  was  by  Stewart  (1971).  If,  in  fact,  the 
phenomenon  reported  by  Watkins  et.  oJL.  was  true  resonance,  an  under- 
standing of  it  would  be  obtained  by  mathematical  analysis  of  simple 
idealized  void  cavity  shapes  via  the  separation  of  variable  method  for 
solving  the  wave  equation. 

Stewart  identified  113  possible  cavity  shapes  describable  through 
coordinate  surfaces  of  the  11  coordinate  systems  in  which  the  scalar 
wave  equation  can  be  separated.  He  ranked  these  surfaces  in  order  of 
mathematical  complexity.  The  simplest  bounded  cavity  shapes  were  found 
to  be  the  sphere,  the  prolate  spheroid,  and  the  oblate  spheroid.  He 
concentrated  his  analysis  upon  these  three. 

The  pertinent  conclusions  of  his  work  were: 

1.  Mathematically,  resonance  cannot  exist  for  void  cavities  peA  Ae. 

2.  If  resonance  of  any  kind  is  to  be  observed  over  a void,  it  has 
to  be  the  interaction  between  the  cavity  and  the  land  surface  above,  not 
a function  of  the  cavity  alone. 

3.  The  resonance  postulated  between  a cavity  and  the  land  surface 
was  not  demonstrated  by  Stewart  to  have  any  physical  meaning  for  spheri- 
cal and  spheroidal  cavities. 

The  conclusion  dr*,wn  from  this  work  was  that  the  development  of  a 
practical  technique  for  the  identification  of  cavities  by  seismic  means 
was  not  likely  to  depend  on  the  measurable  presence  of  resonance,  but 
will  more  likely  be  based  upon  otu^r  wave  phenomena. 


CAVITY  SURFACE  WAVES 


* 


Rayleigh  Waves 

Watkins  OA.  aJL.  proposed  that  the  radial  oscillations  of  a cylinder 
predicted  by  the  work  of  Biot  was  an  adequate  explanation  for  the  reso- 
nance that  they  had  detected  above  a lava  tube.  Since  Blot's  discussion 
of  waves  on  cylindrical  voids  and  the  derivation  completed  for  this 
study  closely  parallel  the  Rayleigh  wave  problem,  it  will  be  useful  to 
look  at  the  assumptions  and  procedures  used  to  solve  this  problem. 

The  question  that  Lord  Rayleigh  formulated  and  answered  was  this: 

If  there  existed  a plane  wave  that  travelled  along  the  plane  interface 
between  an  elastic,  isotropic,  homogeneous  solid  and  a void,  what  would 
be  its  equations  of  motion?  In  order  to  make  this  problem  solvable,  it 
also  had  to  be  assumed  that  the  wave  attenuated  away  from  the  boundary. 
Such  a wave  was  later  found  to  actually  propagate  within  the  earth — the 
Rayleigh  wave.  The  following  discussion  o£  the  Rayleigh  wave  derivation 
is  after  Kolsky  (1963). 

The  boundary  between  the  solid  and  the  void  will  be  the  xy-plane 
with  the  z-axis  positive  towards  the  interior  of  the  solid  (see  Figure 
I).  The  plane  wave  will  move  in  the  positive  x-direction  with  all  part- 
icle motion  restricted  to  the  xz-plane.  By  Helmholtz  Theorem  (Morse 
and  Feshback,  1953)  the  displacement  vector  U can  be  written 

U - V*  + V X ? where  V?  - 0, 
which  becomes  in  two-dimensions 

v*  + <\.  + vi- 

U - U(x*i  + U(t)k, 


U can  also  be  Vt/ittm 


6 


and  this  implies 


U(X)  - ♦ - 4- 

,x  ,z 


U(z)  « ♦ + f 

,z  ,x 


The  wave  equations  in  terms  of  displacement  are 

,(x)  - (X  + p)0  + „v2U(x) 


PU 

PU 


,tt 

(z)  - (X  + p)0  + yV2U 

,tt  ,z 


(z) 


O) 

(2) 

(3) 

(4) 


By  substituting  the  set  of  equations  (1)  and  (2)  into  (3)  and  (4)  respec- 
tively, and  by  grouping  terms,  we  get 


(*  „„  - {(X  + 2y) / p }V24>)  - (T  - (m/p)V2t)  = 0 

{(X  + 2y)/p  }V2$)  - (f  - (y/ p)  V2'f)  - 0. 

, tt  f Z t LL  9x 

These  two  equations  will  be  satisfied  if  and  only  if 


* - {(X  + 2p)/p}V2*  - 0 

, ct 

» tt  ■ <h/p>V2¥  * 0, 
or 

* tt  - C2V2*  where  Cj  - {(X  + 2v)/p}*S  (5) 

» tt  * C2v2t  where  C2  " (6) 
Cj  and  Cj  are  the  velocities  of  the  P-  and  S-wave,  respectively.  The 
general  solutions  to  the  above  equations  are 

♦ » F(z)exp{-t(oit  - kx)}  (7) 

¥ ■ G(z)exp{^(ut  - kx)},  (8) 

when  It  is  assumed  that  ♦ and  ¥ are  harmonic  functions.  It  also  follows 
from  equations  (5)  and  (7)  that 


7 


F(z)  * Aexp{-z(k2  - w2/C 2)^}  + A expfz(k2  - u2/C2)}>4. 

The  term  A exp{z(k2  - u»2/C2)f*  implies  that  F(z),  and  therefore  *,  in- 
creases exponentially  with  depth,  ‘this  is  contrary  to  the  original  as- 
sumption that  the  wave  vanishes  with  depth.  Hence,  to  meet  this  inicial 

* 

condition,  A must  be  zero.  The  formula  for  C(z)  is  abbreviated  in  a 
similar  fashion  and  the  complete  solutions  for  (5)  and  (6)  are 


Aexp{-Qz  + -t(ut 

- kx)} 

where  Q ■ (k2 

- kp*1 

(9) 

Bexp{-Sz  + ,t(u>t 

- kx)} 

where  S * (k2 

- k^. 

(10) 

The  boundary  between  the  solid  and  the  void  must  be  stress-free  and  this 

iplies  that  a.  ^ and  a.  . are  equal  to  zero.  From  the  general  rela- 
Kzz)  (zxj 

tionsbip, 

c(ik)  “ *06ik  + 2Me(ik)’ 


it  can  be  shown  that  in  rectangular  coordinates  that 
o 


(zz) 


(\  + 2 p)$  + A$  + 2 

,zz  ,xx  ,xz 


p(2* 


- y + y 
> * 


.)• 


(ID 

(12) 


(zx)  ' ,xz  ,zz  ,xx' 

Setting  the  left  side  of  these  statements  equal  to  zero  and  substituting 
in  the  designated  partial  derivatives  of  (9)  and  (10),  it  develops  that 


A{  ( A + 2u)Q2  + Ak2}  + B{2pl Os}  0 (13) 

A{2Qlu.}  - B{k2  + S2}  - 0.  (14) 

Substituting  the  second  equation  into  the  first  gives 

4gQSk  - {(A  + 2p)Q2  - Ak2}(k2  + S2).  (15) 


Q - <k2  - k2)** 

S - (k2  - k2)** 


By  recalling  that 


and  expanding,  equation  (15)  is  transformed  to 


C6  - 8r.4  + (24  - I6a 7)r.2  + (16a2  - 15)  = 0 

with 

r.2  - (k2/k)  - (C/C2) 
a2  - (1  - 2a)/ (2  - 2o). 

This  is  a ruble  equation  in  terms  of  r,2,  which  can  be  solved  when  q, 
Poisson's  ratio,  is  known  for  the  solid.  C is  the  velocity  of  the  Ray- 
leigh wave  and  it  is  obtained  from  the  formula  for  r±.  This  velocity  is 
independent  of  frequency,  which  indicates  that  there  is  no  dispersion. 
Hence,  when  values  are  given  for  £,  _A,  and  jj,  it  is  possible  to  calcu- 
late Cj,  C^,  C,  £,  Q,  and  j>.  With  these  quantities  available,  particle 
motion  can  be  studied.  By  inserting  the  designated  partial  derivatives 
of  (9)  and  (1C)  into  equations  (1)  and  (2),  and  noting  that  from  (14) 
that 

B - (2AQia)/(k2  + S2), 

the  fallowing  equations  for  the  components  of  the  displacement  vector 
are  obtained; 

U(x)  » Ak{exp(-Qz)  - 2QSexp(-Sz)/(k2  + S2)  }{sin(u»t  - kx) } (16) 

- Ak-tTexp(-Qz)  - 2QSexp(-Sz)/(k2  + S2)  }fcos((i>t  - kx)} 

l/Z^  * -AQ{exp(-Qr)  - 2k2exp(-Sz)/(k2  f S2)}{cos(u)t  - kx) } (17) 

+ AQ<{exp(-Qz)  - 2k2exp(-Sz)/(k2  + S2)}{ain(ut  - kx) } . 

Both  components  are  complex  numbers  and  the  usual  argument  made  at  this 
point  is  that  since  only  real  solutions  are  desired,  all  values  associ- 
ated with  the  imaginary  term  i can  be  neglected.  This  is  an  invalid 
statement.  As  an  example,  in 
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both  a and  b are  real  numbers,  even  though  this  Is  a complex  expression. 
The  Imaginary  term  J acta  as  a label  that  partitions  the  two  real  num- 
bers a and  b,  maintains  proper  separation  during  mathematical  operations, 
and  controls  their  signs.  In  equat ions ( i6)  and  (17),  both  the  real  solu- 
tions 

U00  ■ Aktexp(-Qz)  - 2QSexp(-Sz)/ (k2  + S2>  HstnOot  - kx)  1 

U(Z)  » -AQ(exp(-Qz)  - 2k2exp(-Sz)/ (k2  + S2)  }{cos(uit  - kx)  I 
and  the  imaginary  solutions 

**  -Ak(exp(-Qz)  - 2QSexp(-Sz)/(k2  + S2)Hcos(u)t  - kx)  } 

U(Z)  ■ AQ{exp(-Qz)  - 2k2exp(-Sz)/(k2  + S2)  }{sin(u)t  - kx) } 
are  correct.  Each  set  of  equations  predicts  the  behavior  of  a surface 
wave  that  has  retrograde  particle  motion  with  respect  to  the  direction 
of  propagation.  The  only  difference  is  that  they  are  90  degrees  out 
of  phase  in  relation  to  each  other.  This  inability  to  discount  the 
imaginary  portion  of  the  equations  at  any  step  of  the  derivation,  except 
perhaps  for  the  last  one,  will  be  an  Important  factor  in  later  deriva- 
tions. 

Axial  Surface  Waves  on  a Cylinder 

Biot  considered  the  possibility  of  a wave  propagating  along  the 
length  of  a bore  located  within  an  elastic,  infinite  medium.  This  dis- 
turbance was  to  be  symmetric  with  respect  to  the  axis  of  the  cylindrical 
void  and  largely  restricted  to  its  surface.  As  with  the  Rayleigh  wave 
problem,  the  particle  motion  was  restricted  to  a plane.  In  this  case, 
it  vas  the  rz-plane  defined  in  the  cylindrical  coordinate  system  (see 
Figure  1).  From  the  Helmholtz  equation  and  the  formulation  of  the  di- 
vergence and  the  curl  in  this  system,  the  following  displacement  vector 
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1b  obtained; 


t 1 


u - (♦  - f )e  + (*  + f/r  + T )e  . 

,r  , z r ,z  ,r  z 

Therefore, 

u(r)  - 4 - <r 

(l) 

,r  ,z 

t/z^  ■ * + f/r  + f . 

,z  ,r 

In  rectangular  coordinates,  the  wave  equation«  are  expressed 

as 

v2*  - (l/c2)*  tt 

(2) 

v2?  - (l/c2)?  tfc. 

(3) 

In  cylindrical  fonn, (2)  becomes 

* + ♦ /r  + ♦ - ♦ /C2. 

,rr  ,r  ,zz  ,tt  1 

Due  to  the  axial  symmetry , the  vector  f is  restricted  to  the 

(A) 

xy-plane 

and  its  components  may  be  written 

(x) 

- -Ysine 

(5) 

- ycose. 

(6) 

The  x-  and  y-components  of  equation  (3)  are 


J*)  + y(x)  + Jx),2  + (x)  M 

,rr  ,r  ,09  ,zz 

Oi 

(7) 

f(y)  + t(y)  + y<y)/r2  + y(y>  _ 

,rr  ,r  ,e9/  ,zz 

(8) 

Substituting  (5)  into  equation  (7)  gives 

f + » /r  - f/r2  + f 

,rr  ,r  ,zz 

’.tt'V 

(9) 

The  same  result  is  obtained  when  (6)  is  substituted  into  (8).  Solutions 
now  oust  be  found  for  equations  (4)  and  (9)  by  separation  of  variables. 


Let 


12 


♦(r,z,t)  • R(r)Z(z)exp(-(wt) . 

The  replacing  of  ♦ In  equation  (4)  with  this  function  results  in  the  fol 
lowing  ordinary  differential  equations; 

7"  - k2Z  - 0 

R"  + R'/r  - (k2  - u>2 / C j ) R - 0. 

The  solution  to  the  first  equation  la  that  of  a wave  travelling  in  the 
positive  z-direction; 

* 

Z • A exp(-Cwz) . 

The  second  differential  equation  has  solutions  in  terms  of  the  Modified 
Bessel  Functions,  Iq  and  (McLachlan,  1955).  Only  approaches  zero 
as  the  radius  goes  to  Infinity.  Since  it  is  assumed  that  the  wave  atten 
cates  away  from  the  boundary,  is  not  considered  in  the  answer.  There 
fore, 

R - A**Kq  ( (k?  - w2/Cj)Jsr) 

and  the  total  solution  to  (4)  is  written 

■ AKQ(mr)exp{<t(kz  - o»t)} 

where 

m ■ k(l  - t;2)*4 

h * v/cr 

The  value  v is  the  phase  velocity.  The  solution  to  (9)  is  found  in  a 
similar  fashion  with  one  of  the  resulting  differential  equations  being 
in  the  form  that  has  solutions  containing  the  Modified  Bessel  Functions, 
1^  and  K^.  At  this  point,  Biot  changed  from  the  exponentials  to  sines 
and  cosines,  making  the  solutions  to  equations  (4)  and  (9) 
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4>  • AKQ(tnr)co8(kz  - wt)  (10) 

V ■ BK  (nr)sin(kz  - u>t)  (11) 

where 


m * k(l 

with  ■ v/C 

n ■ k(i 

- ‘i>’* 

with  t2  * v/C 

Using  these  two  solutions,  properties  of  Modified  Bessel  Functions,  ard 

the  assumptions  that  a,  * and  a,  . vanish  on  the  interface,  it  can  be 

(rr)  (rz) 

shown  that 

(l2) 

4(1  - cpV/(*D/A)U  - < + Kq((ttD/X)  (1  - c2)‘*)/K1((nD/X)(l  - C2)**)  ) 

- 2(2  - C§)  (1  - ^a2)‘S/(nD/X)(l  - t2a2) 

- ((2  - c|)/U  - C2a?)Sa0((irD/A)(l  - Cjt^^j/KJ^D/AMI  - C2^)**)  > - 0, 

where 

a2  - (1  - 2o)/(2  - 2a) 

4 ' v/C2 

D - di&meter  of  cavity 
f ■ wavelength. 

It  also  can  be  shown  tnit  when  the  partial  derivatives  of  equations  (10) 
and  (11)  are  substituted  into  equatiom,  (1),  the  result  is 

U(t)  - A{-mKj  (ma)  - 2k(l  - (aa)/(2  - ;2)  tcos(kz  - «t)  (13) 

U(z)  - A{-kK0(ma)  + 2(l/n  + l)(l  - (ma)/a(2  - t2)  (14) 

+ 2(1  - c2)*1 (ma)K0(na)/(2  - t2 )K1(na}  }sin(k*  - wt), 

where  a is  the  radius  of  the  cavity. 

The  formulas  giver  in  (12),  (13),  and  (14)  arc  important  for  sever- 
al reasons.  They  shoe  the  similarity  between  Blot's  work  and  the  deriva- 
tion for  the  Rayleigh  wave.  They  also  indicate  whether  cylindrical 
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surface  waves  are  applicable  to  the  problem  of  cavity  resonance.  In 
(12),  when  Poisson's  ratio  o is  given,  the  equation  becomes  a function 
of  X/B  and  This  formula  shows  that  unlike  the  Rayleigh  wave,  the 

cylindrical  surface  wave  suffers  dispersion,  with  the  wave  velocity  de- 
pendent on  the  frequency.  As  with  the  Rayleigh  wave,  the  particle  motion 
is  elliptical.  However,  equations  (13)  and  (14)  demonstrate  that  this 
motion  car.  advance  in  relation  to  direction  of  propagation  as  well  as 
be  retrograde.  This  is  a function  of  the  circumference  of  the  cavity. 

Watkins  tt.  at . used  the  equation 

f - C2/1.55D 

to  calculate  the  theoretical  frequency  of  radial  resonance  for  a cavity 
when  dispersion  was  at  a minimum.  This  was  obtained  from  the  plot  of  ^ 
versus  X/D,  whose  values  were  defined  by  (12).  However,  the  components 
of  the  displacement  vector  and  the  original  statement  of  the  problem  both 
indicate  that  this  is  not  a standing  wave  phenomenon,  but  a wave  that 
stoves  parallel  to  the  z-axis.  As  with  the  Rayleigh  case,  to  have  waves 
possessing  set  frequencies  sweep  past  a single  point  over  an  extended 
period  of  time,  it  is  necessary  for  these  waves  to  be  continuously  gener- 
ated. This  is  not  the  case  for  energy  dispersed  by  a dynamite  explosion. 
Consequently,  Biot's  derivation  cannot  explain  the  definite,  long-lived 
oscillations  recorded  by  Watkins  tt.  at.  above  some  cavities.  A different 
interpretation  must  be  made. 

Circumferential  Surface  Waves  on  a Cylinder 

Since  resonance  was  found  over  relatively  shallow  caves  and  because 
surface  waves  attenuate  away  from  the  wall  of  a cavity,  it  was  thought 
that  a different  type  of  surface  wave  might  be  the  solution  to  the  problem. 
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It  was  believed  that  resonance  may  be  Induced  bv  cylindrical  surface 
waves  travelling  about  the  circumference  of  the  cavity.  With  this  form 
of  wave,  the  frequency  and  strength  of  oscillations  recorded  at  the 
earth's  surface  would  be  a function  uf ; t he  wave's  velocity,  the  cavity's 
circumference,  the  depth  to  the  cave,  and  the  types  of  hit erl erenoe  pro- 
duced by  the  different  surface  waves.  If  this  theory  were  correct.  It 
would  mean  that  a great  deal  could  he  told  about  the  cavity  by  the  traces 
of  the  seismogram.  It  would  also  explain  why  resonance  had  not  been  re- 
corded over  deep  cavitieR. 

The  statement  of  the  problem  Is  similar  to  the  previous  two  exam- 
ples. Tt  is  assumed  that  there  is  a cylindrical  bore  in  an  elastic.  In- 
finite solid,  which  has  a symmetric  wave  moving  about  its  circumference. 
The  particle  motion  is  restricted  to  the  r9-plane  and  it  attenuates  away 
from  the  surface  of  the  cavity. 

By  writing  the  divergence  and  the  curl  of  the  Helmholtz  Theorem  in 
terms  of  cylindrical  coordinates  and  invoking  the  condition  that  no  mo- 
tion can  occur  outside  of  the  rO-plune,  the  components  of  the  displace- 
ment vector  U are  expressed  as 

U(r)  - ♦ + y /r  (1) 

• r 

U(9)  - ♦ /r  + T . 

,6  ,r 

After  Eringer  (1967),  the  most  general  expression  of  the  wave  equation  is 

pU  ..  - (X  + 2p)V0  - u(?  X V X U).  (2) 


The  dilatation  0 is  given  by 

G - V’U 


U(r)/r  + U<r)  V(V/r 
* r » 
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0 - (*>r  ♦ »yr )/r  + (*>rr  - V’  + »f0r/r> 

* <*,99/r  - \or,'r 


n ■ /r  + ♦ + <f  /r* 

u \r  ,rr  ,oo 


0 3 v2*. 

This  demonstrates  that  the  dilatation  Is  equal  to  the  Lap lac  Ian  of  the 
scalar  potential  $.  The  curl  of  II  is  given  by 


v X V X U - (U(0)/r2  + U(0Vr  - U^/r^e,  + (U(tt)/r2  - U(0)/r  (4) 

•e  » 00  1 * r 

- u(0)  - 11^/r2  + U^/r)e  . 

,rr  ,6  ,re  0 

By  substituting  equation  (3)»  equation  (4),  and  the  partial  derivatives 
of  the  components  of  the  displacement  vector  into  (2),  the  following 
equations  are  obtained; 


pU(*J  « (X  + 2p)(V2*)  + (p/r)  (f  + » Jr  + t / r2) 

,tt  ,r  »rr  ,r  ,uo  ,u 

pU(0J  - {(X  + 2p)/r) (V24>)  0 -»*(»  + » h + » ee/r2)  . 

By  inserting  the  values  of  and  U^0^  given  in  (1)  and  observing  that 

the  Laplecian  of  t is  present  on  the  right  side  of  the  above  equations, 
we  get 

p(*  + (p/r) (1»  ) - (X  + 2p)(V2*)  + (p/r) (v2y) 

,tt  ,r  ,tt  ,u  ,r  ,o 

(p/r)  (♦  ) A - p(f  .J  r - ((X  + 2p)/r}(V2*)  - p(72f) 

itt  ,0  #tc  tr  »0 

or 

(*  tt  - {(X  + 2p)/p}V2*)>r  + (l/r)(f>tt  - (p/p)72f)0  - 0 
d/r)(«  - {(X  + 2p)/p}V2*)  - (t  - (p/p)72f)  - 0. 

|CC  |0  »4 

These  will  be  valid  if  and  only  if  the  following  equations  hold; 

♦ tt  “ c2?2*  where  Cj  - {(A  + 2p)/p}>S 


(5) 
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i. 

y tt  " where  C2  - (u/p)  (6) 

These  sre  the  two  scalar  wave  equations.  We  have  shown  that  they  are  the 
same  in  cylindrical  coordinates  as  these  derived  in  the  Rayleigh  wave 
problem.  The  next  step  is  to  solve  equation  (5)  by  separation  of  vari- 
ables. The  answer  to  equation  (6)  will  be  of  the  same  form.  Let  ♦ have 
a solution  that  is  harmonic  with  time.  Hence, 

♦ ■ R(r)0(Q)exp(-tut ) . 

Substituting  into  (5)  results  in  the  expression 

V2(RO)  ♦ k2R0  - 0 where  k - u)/^. 

Expanding  the  Laplacian  in  cylindrical  coordinates  gives  the  equation 

r2R  /R  + rR  /R  + r2k2  - -0  .JO  - m2 
,rr  tr  i j u □ 

where  m is  the  integration  constant.  This  can  be  separated  into  two  or- 
dinary differential  equations; 

e,99  + «2e  ‘ 0 

r2R  + rR  + (r2k2  - m2)R  - 0. 

,rr  ,r  1 

The  first  equation  has  the  general  solution 

* * 

0 ■ A exp(-xjnS)  + B exp(<un6). 

The  negative  exponent  term  is  for  a wave  moving  about  the  cavity  counter- 
clockwise and  the  positive  exponent  is  for  one  moving  clockwise.  Choosing 
the  counterclockwise  wave  gives  the  solution 

★ 

0 « A exp(-tm0) . 

The  second  equation  is  in  the  form  of 

z2w,zz  + zw,z  + (z2  " y2)w  " °» 

which  is  the  differential  equation  for  the  Hankel  Function.  Therefore, 
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the  general  solution  for  the  second  differential  equation  is 

**  (2)  **  (H 

R = A IT  ;(rk.)  + B ||v  ;(rk,). 
ml  ml 

According  to  Butkov  (1968),  when  exp(tait)  is  chosen  as  the  lime  harmonic 
portion  of  the  original  function,  it  is  necessary  to  choose  the  tiankel 
Function  of  the  second  kind  so  as  to  satisty  the  Sommer  fold  Radiation  con- 
dition. Since  this  is  the  only  Hankel  Function  used  throughout  the 
derivation,  the  superscript  will  be  left  off.  The  solution  for  R is 


now  written 


irk 

R * A H (rk,) 
m 1 


with  the  Hankel  Function  of  the  second  kind  defined  as 


H (rk.)  = J (rk.)  - iY  (rk.) 
ml  ml  ml 


Therefore,  the  solutions  for  and  T are 


<J>  = AH  (rk.  )exp{-t(u>t  - m0)} 
m i 


¥ = BH  (rk_)exp{-t.(iot  - m6)}, 
m i 


where 


kt  = u/C 

k 2 = w/C^. 

Expressing  the  stress  components  °(rr)  and  °(re)’  *n  terms  ♦ and 
^ results  in  the  following  two  expressions;  ' 

°(rr)  " U + 2“>(*,rr  " \eh2  + ',8r/r)  + (,) 

"\Jr2  - \re/r  + \r'r  + *,8/r2> 

"(re)  • <u/r2)(2,‘,r0  + ’.89  + r*,r  ' r2’,rr  ' 2»,9>’  (l0) 

The  first  derivative  of  the  Hankel  Function  is 

V'V.r"  <V2>(Vl<rkl>  - Vl(rkl»- 
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The  second  derivative  is 


H (rk  ) - (k. / A ) { H (rk,)  - 2H  (rk  ) + H (rk  )}. 

m l ,rr  1 m— £ I m l nrrz  ] 

Using  these  derivatives,  it  is  now  possible  to  take  the  part  inis  of  (7) 

and  (8)  and  substitute  them  into  equations  (9)  and  (10).  For  simplicity, 

we  will  adopt  the  following  convention; 


- H (rk, 

m 1 

- H (rk  ) 

III  4 


The  stress  components  now  become 


°<rr>  ■ (*«>*k5/2  + xkl/4)(\-2  ' 2H„  + "*2>  - <x"2/r2)H» 

(ID 

+ (Xk,/2r)(H  , - H , ) } + B{(2yroe/r2)h 

i m- 1 nrr  l m 

t 

- (pnw.k2/r)(hm_1  - h^j)  })expU(wt  - m6)} 

*v 

°(re>  ■ ("«->‘k2/4>(V2  - 2hm  + W ‘ 

(12) 

+ (pk2/2r)(h^1  - h^j)}  + A{  (2pnvt/r2)Hm 

- (umtkj /r)  (Hml  - H^j)  ))exp{-t(u>t  - me)}. 

If  «i  is  the  radius  of  the  cavity,  the  stresses  expressed  above  vanish  when 
r - a.  Given  radius  a and  letting 


o ■ (yk2/2  + Ak2/4)(H  „ - 2H  + H ) - (Am2/a2)H 

l l D'<£  m nri-z  m 

+ - Vi> 

6 - <k|/4)  (ha_2  - 2b.  + - (k^aXh..,  - h^,) 

y - <2W.*W.  - - H^,) 

6 » (2ym/a2)hm-  (ymkj/a)  (h^  - h^), 


equations  (11)  and  (12)  become 


r 
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0 ■ Ao  + B6i 
0 ■ BB  + A y-i. 

Substituting  the  second  equation  into  the  first  results  in 

u£  ■ 

Let 

P - (yk2/2  + Xk2/i)(J  - 2J  + J ) - (Xra2/a2)J 
l l m-  / m nn*t  n 


(13) 


m 


+ (Xkl/2a)(Vl  - *W 


Q - (yk2/2  + Ak2/4) (Y  0 - 2Y  + Y ,)  - (Am2/a2)Y 
l i m-/  m niT4  m 


+ '“■““''Vi " Vi> 

R - <k!'*><V2  - 2^  + W + ("2/a2>J„  - <V2a><Vl  - W 

S - <^M><V2  - 2jm  + V2)  + (m2/a2)y|n  - <V2a)(Vl  * Vl> 

T - (2n,/a’)Jm  - (-^/a)  (J^,  - J^,) 

U - (2m/a2)Y  - (mk  /a)(Ym  , - Y . ) 

in  l m— i m+i 

V - (2yo/a2)jm  - (ymk2/a)  (j^  - j^) 

W - (2ym/a2)y  - (ymk  /a)(y  - 

in  / m-i  pit  i 

Expanding  o,  Ji,  and  A into  their  real  and  imaginary  parts  creates  th. 
following  complex  formulae; 

a " P - Q-t 
6 - R - Si 
Y - T - Ui 
6 - V - Wi. 

Substituting  these  into  equation  (13)  gives 

(PR  - QS)  - (QR  + PS)i  - (TV  - UV)  - (UV  + TW)i. 


These  two  complex  expressions  will  be  equal  to  one  another  if  and  only 
If 
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PR  - QS  » TV  - UW 


(H) 


and 


QR  + PS  * UV  + TW.  (15) 

Therefore,  the  u>  that  causes  both  of  these  equations  to  hold  simultan- 
eously, will  be  the  angular  frequency  of  the  wave  propagating  about  the 
circumference  of  the  cavity.  Using  this  frequency,  it  is  possible  to 
describe  the  particle  motion.  Inserting  the  partial  derivatives  of 
equations  (7)  and  (8)  into  equations  (l)  gives 

U(r)  - {(Ak./2)(H  , - H . ) - (Bm-t/r)h  }exPU(ujt  - me)}  (16) 

l m- 1 m+ 1 m 

U<0)  - - (Bk2/2)(hm_1  - Vl*  }expU(wt  - me)}.  (17) 

We  already  have  that 

0 ■ B8  + Ay*. 


or 

B - A(y/B)^, 

which  Indicates  that 


B - A({(UR  - ST)/(R2  + S2)}  + {(TR  + US)/(R2  + S2)U). 

Replace  B in  equations  (16)  and  (17)  with  the  above  value  and  expand  the 
Hankel  Function  into  its  real  and  imaginary  parts.  This  changes  the  com- 
ponents of  the  displacement  vector  to 

U(r)  - A({rcos(ut' - mB)  - Asin(ut  - m0)} 

-**s 

+ (rsin(o>t  - m6)  + Acos(u>t  - m8>}<) 

(e ) t 

U ■ A((5cos(u)t  - me)  - nsin(ut  - m0)} 

+ {Esin  (out  - m6)  Qcos(o)t  - m0 )}.(.) 


where 
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(k./2)<1  . - j - (m/ a) { (UR  - ST)/(R2  + S 2)}y 

I m-1  mrl  m 

+ (m/a) { (TR  + US)/(R2  + S2) }j 

m 

(-k^KY^  - Y^)  - (m/a)  { (UR  - ST)/(R2  -f  S2)}.jm 

- (m/a) { (TR  + US)/(R2  + S2))j 

m 

<~m/a)Y  - (k  /2){(UR  - ST)/(R2  + S?) } (J  . - J ) 

m i m-i  nrri 

- (k2/2){(TR  + US) / (R2  + S2)}(ym  l - y^j) 

(-m/ a) Jm  + (k2/?.){(’JR  - ST)/ (R2  + S2)}(ym  l - y^) 

- (k2/2){(TR  + US)/(R2  + S2)}(jtn  l - j^). 

S,  T,  and  U are  as  defined  earlier  and  a is  the  radius  of 
It  is  not  possible  to  find  the  absolute  values  for  the 
displacement  vector  U,  but  since  both  components  are  expressed  in  terms 
of  A,  their  amplitudes  may  be  compared.  If  the  constant  A is  negative  or 
complex,  the  only  change  to  and  would  be  a phase  change.  There 

also  is  no  loss  of  generality  if  only  the  real  portion  of  the  above  equa- 
tions is  used.  Therefore,  the  components  of  the  displacement  vector  may 
be  witten 

- A{rcos(uit  - m0)  - Asin(u)t  - m0)}  (18) 

(8 1 

Uv  * A{Ecos(ut  - m0)  - ftsin(ut  - n»8)}.  (19) 

This  completes  the  derivation  of  the  circumferential  surface  wave  on  a 
cylinder.  The  Intention  was  to  use  these  displacement  formulae  to  define 
the  motion  of  the  particles  along  hypothetical  horizons  above  the  cavity 
as  showu  in  Figure  2. 

A computer  program  was  written  to  numerically  solve  equations  (14) 
and  (IS),  when  a,  X.,  and  ^ were  given.  The  procedure  was  to  set  m,  the 
mode,  end  vary  w iteratively  until  the  left  side  of  the  equations 


r - 
a * 


a 

The  values  R, 
the  cylinder. 
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Figure  2. 


Possible  elllpticsl  action  end  direction  of  rotation  of  points 
along  two  hypothetical  lines  above  a cavity. 
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equalled  the  right.  It  was  found  to  be  impossible  to  satisfy  both  equa- 
tions simultaneously.  Even  when  imaginary  equation  (15)  was  neglected, 
the  angular  frequency  generated  for  a cavity  of  reasonable  size, resulted 
in  a wave  whose  velocity  was  much  greater  than  that  of  a P-wave.  Since 
it  was  assumed  in  the  derivation  that  the  cavity  surface  wave  would  be 
a combination  of  compressional  and  shear  motion,  it  would  he  impossible 
for  a wave  to  propagate  at  this  speed. 

This  conclusion  is  in  agreement  with  St  like,  who  completed  a three- 
dimensional  version  of  Blot's  problem.  He  found  that  the  cylindrical 
surface  wave  could  spiral  about  a tube  as  illustrated  in  Figure  3,  but 
that  for  all  modes  except  m » 1,  no  wave  could  have  only  r6-motion.  For 
m « 1,  this  occurred  when  the  wavelength  approached  infinity  or  the  cyl- 
inders radius  approached  zero.  The  results  of  our  derivation  and  those 
of  Stilke  imply  that  circumferential  surface  waves  on  a cylinder  cannot 
exist  and  therefore  cannot  be  the  cause  of  apparent  cavity  resonance. 


FLEXURAL  WAVES 


With  the  failure  of  cavity  surface  waves  to  explain  cavity  reso- 
nance, the  next  hypothesis  investigated  was  the  "drumhead"  effect  as  sug- 
gested by  Finn  Bronner  (personal  communication).  The  existence  of  such 
an  effect  would  imply  that  the  roof  of  a cave  oscillates  vertically  with 
little  or  no  horizontal  components  of  motion.  As  a first  step  In  studying 
this  possible  phenomenon,  the  behavior  of  a vibrating  membrane  was  exam- 
ined. 


Vibrating  Membrane 

The  equation  of  motion  for  a round  membrane  in  cylindrical  coor- 
dinates is 


where 


i ,r  + 5 Jr  - rH  - 0 

,rr  ,r  ,tt 


(1) 


r2  - m/T, 


with 


m ■ mass  per  unit  area 
T ■ total  tension. 

Using  the  boundary  condition  for  £(r,t)  of 

C(a,t)  - 0 

for  all  time  Jt  and  membrane  radius  a,  and  the  lntitial  condition 

5(0,0)  - eQ. 

it  can  be  shown  that 


J0(wa/D  - 0. 


The  radial  frequency  of  the  membrane  will  be  that  w which  satisfies  this 
equation  (Churchill,  1963).  There  are  a number  of  implicit  assumptions 
that  have  to  be  made  to  solve  this  proble.  They  are; 
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1)  vibrations  occur  within  a vacuum, 

2)  there  is  no  damping, 

3)  the  system  is  perfectly  elastic, 

4)  amplitudes  are  small, 

5)  deformation  due  to  gravity  may  be  neglected, 

6)  the  membrane  is  thin  (surface  of  "points"),  and 

7)  tension  is  constant  at  all  times,  even  during  deformation. 

The  first  five  assumptions  present  no  problems,  since  they  are  used  ex- 
tensively in  wave  propagation  problems  whose  results  have  been  shown  to 
be  reliable.  The  last  two  assumptions  cannot  be  met.  No  cavity  roof  can 
be  thought  of  as  being  a surface  of  points  and  the  forces  that  support  a 
roof  are  not  analogous  to  tension.  Therefore,  equation  (1)  is  not  ap- 
plicable for  calculating  frequencies  found  over  caves.  The  next  approach 
was  to  look  at  the  equations  or  a vibrating,  thin  plate. 

Vibrating  Plate 

The  general  equation  of  motion  for  a thin  plate  is 

v4e  + <1/C2K  - 0 (2) 

» c t 

where 

C - (Eh2)/I2p(l  - a2),  (3) 

with 

E *>  Young's  Modulus 
h ■ thickness  of  plate 
p ■ density 
o » Poisson's  Ratio. 

In  the  solutions  that  follow,  all  the  assumptions  are  the  same  as  for  the 
membrane  problem,  except  that  the  last  two  are  replaced  by  one  assumption 
— the  plate  is  thin.  What  this  means  is  that  within  the  plate  there  is 
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a neutral  surface  which  is  neither  compressed  nor  extended  as  the  plate 
bends,  and  that  the  motion  of  this  hypothetical  surface  closely  mirrors 
the  movement  of  the  plate's  top  and  bottom  surfaces.  There  is  no  exact 
definition  of  when  a plate  can  no  longer  be  considered  thin,  but  perhaps 
it  is  when  the  plate's  behavior  can  not  be  predicted  theoretically.  In 
seismology,  the  oscillations  associated  with  plates  are  called  flexural 
waves,  and  they  are  found  to  exist  over  frozen  lakes  (Ewing,  Jardctzky, 
and  Press,  1957)  and  in  desert  environments  where  a crust  of  caliche  has 
formed  (Richter,  1958).  The  question  is  whether  the  resonance  found  by 
Watkins  it.  at.  over  a lava  tunnel  and  at  an  underground  nuclear  test 
site  was  due  to  the  generation  of  flexural  waves. 

To  solve  equation  (2)  for  a circular  plate,  the  boundary  conditions 
for  £(r,0)  are  made  to  be 


S(a,8)  - 0 

Ua,0)  r - 0, 

9 r 

where  a is  the  radius  cf  the  plate.  The  result  is 


which  has  n solutions  for  each  Integer  m.  X is  found  by  a series  ex- 
pension  and  is  constant  for  set  values  of  m and  n.  The  equation  for  the 
angular  frequency  is 

» „ - (hX  /a2){E/\2p(l  - o2)}**.  (4) 

m,n  mtn 

Figure  4 has  some  examples  of  modes  of  vibration  for  a circular  plate 
(McLachlan,  1951). 

Landau  and  Llfshitz  (1959)  did  a similar  problem  for  a rectangular 
plate  with  sides  of  length  a and  b.  The  boundary  conditions  are 

C(x.y)  - 0 

for  x ■ 0 and  a 

e(x,y),x  ■ 0 


Figure  a.  Nodes  of  vibretlon  for  clrculer  end  rectangular  pletes.  Patterns 
bacons  progreaaively  no re  coeplea  ea  the  mode  increases.  Higher 
nodes  ere  rarely  seen  because  they  attenuate  very  quickly. 
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£(*.y)  - o 

C(x,y)  - 0 

• 7 


for  y ■ 0 and  b. 


The  solution  to  equation  (2)  is 

(,  - Asin(8*x/a)sin(txy/b) , 

where  s and  t are  integers.  The  frequencies  are  given  by 

f - a)  /2n 
s,t  8,t 

f . - (hir/2)(s2/a2  + t2/b2){E/12p(l  - o2)}\  (5) 

8,t 

Since  the  tunnel  beneath  the  Pisgah  lava  flow  had  a ceiling  that  could  be 
approximated  by  a rectangle,  the  characteristic  vibrations  were  found 
using  quatlon  (5).  The  following  were  the  values  needed; 

Poisson's  Ratio  o ■ 0.25 

toung's  Modulus  E ■ 2.33  X 109  nt/m2 

Density  p ■ 2.20  X 103  kg/m3 

Length  a » 100.0  m 

Width  b *•  10.0  m 

Thickness  h * 2.5  m. 

The  constants  a,  a,  b,  and  h were  taken  from  Watkins  et.  at, . A stan- 
dard value  for  vesicular  basalt  was  assigned  to  E (Clark,  1966).  Sub- 
stituting all  the  above  values  into  equation  (5),  and  setting  a and  t_ 
equal  to  one,  gives 

f.  - (2. 5w/2.0)  (1/101*  + 1/102)  {2. 33  X 109/12(2.20  X 103)  (1  - .252)}1* 
fl  j • <3.91)(1.01  X 10”2) (4.89  X 102) 


1.1 


19.3  hertz. 


Higher  modes  result  in  higher  frequencies  and  as  shown  In  Figure  4,  more 
complicated  patterns.  Higher  modes  are  often  not  found,  since  they  tend 


to  dampen  quickly.  The  above  answer  is  the  simplest  characteristic  vi- 
bration for  a rectangle  and  it  is  very  close  to  the  16  hertz  signal  as- 
sociated with  the  Pisgah  lava  tunnel. 

The  information  on  the  dimensions  of  the  chimney  formed  by  an  un- 
derground nuclear  explosion  and  on  the  rock  type  capping  it  was  not  exact 
enough  to  satisfactorily  work  the  frequency  problem.  However,  by  making 
several  different  sets  of  assumptions  for  equation  (4),  the  frequencies 
calculated  usually  came  to  be  between  10  and  20  hertz.  The  field  mea- 
surement was  12  hertz.  The  data  for  the  Kana-a  lava  tunnel  were  insu- 
ficient  to  even  make  this  kind  of  approximation. 

In  the  previous  section,  we  proved  that  neither  axial  or  circumfer- 
ential cavity  waves  could  explain  the  resonance  found  over  a few  under- 
ground chambers.  In  this  section,  we  have  shown  that  the  frequencies  of 
flexural  waves  agree  well  with  those  measured  in  the  field.  This  Implies 
that  flexural  waves  may  be  of  value  in  detecting  certain  types  of  caves. 
The  special  characteristic  needed  for  their  generation  is  a thin  roof 
with  a large  surface  area.  This  feature  would  explain  why  resonance  has 
not  been  found  by  Godson  and  Watkins,  Rechtien  and  Stewart,  and  Stllke 
for  relatively  deep  caves.  Presently  this  ic  an  unproven,  but  promising 
theory.  To  conclusively  prove  that  cavity  resonance  is  due  to  flexural 
waves,  it  will  take  extensive  field  work  in  comparable  areas,  with  spe- 
cial attention  to  cavity  dimensions  and  the  physical  properties  of  the 
surrounding  rock. 


1^ 

RAY  TRACING 

Cook  found  a decrease  in  the  normalized  amplitudes  of  waves  re- 
corded above  brine  cavities.  He  concluded  that  these  anomalies  were 
caused  by  a seismic  "shadow"  formed  when  the  cavity  acted  as  an  obstacle 
to  a portion  of  the  up-coming  reflections.  Watkins  zt.  aJt.  had  a similar 
amplitude  decrease  above  the  Kana-a  lava  tunnel  and  over  the  U4B  under- 
ground nuclear  test  Jite.  At  this  location  they  also  recorded  delays  in 
arrival  times.  Apparent  delay  times  were  also  seen  over  some  Missouri 
caves  investigated  during  thi<:  study.  The  first  attempt  to  theoretically 
study  this  "shadow"  phenomenon  was  with  ray  tracing. 

Cavanaugh  (1973)  had  developed  an  efficient  computer  procedure  for 
the  iterative  tracing  of  rays  through  volumes  containing  complicated 
boundaries.  This  differed  from  the  usual  methods  by  being  able  to  handle 
high-powered  polynomial  surfaces  as  easily  as  simple  surfaces  and  by 
using  the  three-dimensional  form  of  Snell's  Law.  This  is  a vector  equa- 
tion and  is  written 

(1/V0)(I  X N)  - (1/VjXR  X N),  (1) 

where 

Vq  “ velocity  of  the  medium  containing  the  incident  ray, 

- velocity  of  the  medium  containing  the  refracted  ray, 

I - unit  incident  ray  vector, 

R - unit  refracted  ray  vector,  and 

N ■ unit  vector  normal  to  the  interface  between  the  two  media. 

Figure  5 shows  a generalized  boundary  intersected  by  the  incident  vector 
X.  The  plane  of  refraction  is  defined  by  I X N.  When  the  magnitudes 
of  the  cross-products  in  (1)  are  calculated,  the  result  is  the  familiar 
expression  of  Snell's  Law  in  two-dimensions. 
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Using  this  vector  equation,  Stravroudla  (1972)  derived  the  fol- 
lowing formula  for  the  unit  refracted  vector; 

R - TN  + p I 

where 

» - Vv0 

T ■ -ycos(i)  + {1  - vi2sin2(i)  J*5. 

When  reflection  Is  considered  as  a special  case  of  reflection,  it  Is  seen 
that  p*l  and  r*-cos(i).  This  implies 

r ■ -2cos(i)N  + I, 

where  r is  the  unit  reflected  vector.  By  defining  the  plane  of  refraction 
and  reflection  at  every  point  of  Intersection  and  by  using  the  equations 
for  R and  r.  It  is  possible  to  mathematically  construct  a ray's  path 
through  almost  any  volume. 

Most  natural  cavities  are  smaller  than  the  wavelengths  of  waves 
that  propagate  through  the  earth's  crust.  Ray  tracing  cannot  be  used  for 
this  case  because  of  dispersion.  However,  it  was  hoped  that  if  could  be 
applicable  to  underground  chambers  with  extremely  large  dimensions,  such 
as  some  brine  cavities. 

In  order  to  produce  the  most  realistic  results,  modifications  were 
made  to  the  basic  ray  tracing  program.  Boundary  phase  changes  were  added 
and  amplitudes  were  adjusted  at  each  point  of  reflection  and  refraction 
so  as  to  agree  with  the  division  of  energy  predicted  by  the  Zospprits 
equations.  However,  this  approach  was  soon  abandoned.  The  exponential 
increase  in  rays  due  to  phase  conversions  for  both  refractions  and  re- 
flections created  tins  consuming  computer  runs,  even  for  the  simplest 
modals.  Ray  tracing  alao  could  not  predict  interference  patterns,  dif- 
fractions, or  the  form  of  any  resulting  surface  waves.  What  was  being 
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done  was  an  attempt  to  take  some  of  the  simplified  properties  of  waves 
(Snell '8  Law,  Zoepprltz  equations,  and  phase  changes)  and  use  them  to 
model  in  detail,  a very  complicated  form  of  wave  propagation.  It  became 
clear  that  this  type  of  modelling  could  only  be  done  successfully  if  you 
went  back  to  the  source  of  these  properties— 'the  wave  equation. 


WAVE  SCATTER 


Because  most  caves  are  smaller  Chan  the  lengths  of  seiomic  waves. 

It  is  not  possible  for  them  to  be  opaque,  even  with  respect  to  a pure 
shear  wave.  For  a large  cavity,  energy  will  be  transferred  to  the  lee- 
ward side  by  various  surface  waves  and  by  diffraction.  Therefore,  in- 
stead of  a cave  having  a seismic  "shadow",  it  will  probably  have  an  area 
of  complex  "turbulence",  which  could  be  compared  to  that  found  behind  a 
pier  pile  after  an  ocean  wave  has  passed.  This  "turbulence"  is  formed  by 
a combination  of  reflected,  refracted,  and  diffracted  waves,  and  it  is 
called  wave  scatter.  To  solve  scatter  problems,  it  is  necessary  to  use 
the  wave  equation,  that  can  be  solved  by  either  analytical  or  numerical 
methods. 

Analytical  Solutions 

Ying  and  Truel  (1955)  worked  the  problem  of  a plane  P-tiave  moving 
in  an  infinite,  elastic  solid,  and  incident  upon  a spherical  obstacle. 
They  considered  three  types  of  spheres;  void,  elastic,  and  rigid.  For 
the  last  two  obstacles,  three  separate  waves  were  derived.  They  were  the 
unimpeded  incident  wave,  the  scatter  wave  formed  at  the  obstacles  sur- 
face, and  the  wave  generated  within  the  sphere.  The  void  would  be  anal- 
ogous to  an  air-filled  cave,  which  unlike  the  rigid  and  elastic  spheres, 
would  not  form  on  internal  wave. 

Ying  and  Truell  completed  all  equations  of  motion  and  set  all 
boundary  conditions  in  spherical  coordinates.  They  also  restricted  part- 
icle motion  to  the  two-dimensional  p0-plane.  The  scalar  potentials 
for  the  three  separate  obstacles  became  expressed  in  Spherical  Bessel 
Functions,  and  thesa  vara  usad  to  calculate  the  scattering  cross-sections 
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for  each  case.  The  components  of  the  displacement  vector  were  not  found, 
but  they  were  obtainable  from  the  equations  provided. 

White  (195^  did  a similar  problem  for  a cylindrical  scatterer  in- 
bedded  in  an  infinite,  elastic  medium.  The  obstacles  studied  were  elastic, 
void,  and  fluid-filled  cylinders.  The  second  and  third  could  be  used  as 
models  for  air-filled  caves  and  ones  containing  water  or  mud.  The  inci- 
dent plane  P-  and  S-wave  possessed  a single  frequency  and  could  strike 
the  cylinder  at  any  angle  of  incidence.  The  vector  potential  was  not  con- 
verted to  a scalar  by  limiting  particle  motion  to  a single  plane.  The 
problem  was  solved  in  three-dimensions,  with  the  scalar  and  vector  poten- 
tials being  expressed  in  terms  of  Bessel  and  Hankel  Functions.  This  der- 
ivation also  demonstrated  that  an  incident  P-  or  S-wave  will  result  in 
scattered  waves  of  both  types.  This  is  referred  to  as  mode  conversion. 

The  derivations  by  White  and  by  Ying  and  Truell  provide  much  infor- 
mation concerning  the  problem  of  scatter,  but  both  have  serious  limita- 
tions, which  restrict  their  usefulness  for  the  theoretical  study  of  cave 
detection.  In  seismology,  ground  motion  is  measured  at  the  earth's  sur- 
face, which  indicates  that  for  the  cave  problem  there  are  two  important 
boundaries;  the  surface  of  the  cavity  and  the  surface  of  the  ground. 

Neither  of  the  above  derivations  considered  a second  planar  boundary.  To 
have  included  such  an  interface,  would  have  created  extreme  difficulties 
with  boundary  conditions,  since  a plane  that  does  not  go  through  the  or- 
igin is  not  easily  expressed  in  cylindrical  or  spherical  coordinates. 
Without  a plane  boundary,  the  interaction  between  the  scattered  wave  and 
the  surface  of  the  ground  cannot  be  accurately  known.  The  second  drawback 
la  that  the  cavities  must  mirror  the  coordinate  system  used.  This  results 
in  targets  of  very  simple  shapes.  Again,  this  is  dictated  by  the  fact 


that  uncomplicated  boundary  conditions  are  needed  to  solve  the  analytical 
problem  efficiently  and  accurately.  These  difficulties  can  be  avoided  by 
solving  the  scatter  problem  numerically. 

Numerical  Solutions 

Both  the  finite  element  and  the  finite  difference  methods  may  be 
employed  to  describe  the  behavior  of  scattered  waves.  These  two  numeri- 
cal procedures  have  some  common  properties.  Except  for  trivial  cases, 
each  requires  the  use  of  computers.  This  has  restricted  their  applica- 
tion in  the  past.  Both  techniques  represent  the  geological  profile  with 
a set  of  points  and  both  use  the  basic  form  of  the  wave  equation; 

pU  - (X  + m)V0  + pV2U. 

f tt 

These  properties  allow  the  problem  to  be  solved  in  the  rectangular  coor- 
dinate system  and  they  permit  very  irregular  cavities  without  difficul- 
ties with  boundary  conditions. 

With  the  finite  element  method,  each  point  represents  a "lumped" 
mass,  and  every  pair  of  points  is  assigned  elastic  properties  that  are 
defined  by  the  types  of  rocks  encountered  and  the  boundary  locations.  To 
generate  wave  motion  through  this  maze  of  mathematical  weights  and  springs, 
the  second  ordered  differential  wave  equation  is  stated  in  terms  of  Hook's 
Law.  A large  group  of  simultaneous  equations  is  formed  when  the  values 
at  the  points  are  substituted  into  the  modified  wave  formula.  These  can 
be  solved  at  set  Increments  of  time  by  the  Runge-Kutte  algorithm.  Each 
set  of  answers  gives  the  positions  of  the  model  points  in  that  time  frame 
and  the  sum  of  these  positions  over  a length  of  time  defines  the  wave 
motion. 

When  using  the  finite  difference  procedure,  the  points  of  the  model 
must  be  spaced  at  regular  intervals.  The  values  assigned  to  each  point 
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depend  on  the  density  and  the  elastic  moduli  of  the  material  surrounding 
that  point.  The  only  change  in  the  wave  equation  is  the  replacement  of 
the  differentials  by  difference  approximations.  These  formulae  are  in  a 
form  that  can  be  easily  programmed.  The  motion  is  found  by  iteration, 
where  the  movement  of  a point  is  directly  dependent  on  the  movement  of 
that  point  and  the  surrounding  points  in  the  previous  time  increments.  I 
have  chosen  to  work  with  this  method  because  it  is  relatively  easy  to  pro- 
gram and  it  is  extensively  discussed  in  the  literature. 

There  are  several  finite  difference  methods  used  to  define  wave  mo- 
tion. The  first  one  to  be  applied  to  the  problem  of  scatter  was  devel- 
oped by  Claerbout  (1970).  He  devised  a good  one-dimensional  approximation 
to  the  two  dimensional  wave  problem.  This  is  a very  attractive  feature  in 
that  it  saves  large  amounts  of  computer  storage  space.  This  is  usually 
the  biggest  problem  encountered  with  finite  difference  methods. 

Claerbout  reformulated  the  wave  equation  in  terms  of  pressure  so 

that 

V2P  - (i/c2)Ptt  (1) 

where' 

C - (x/p)1* 
k - \ + (2/3)p. 

The  term  jc  stands  for  incompressibility.  The  behavior  of  this  pressure 
wave  is  idrntfoii  to  t,  , of  the  P-wave.  The  solution  to  equation  (1) 
in  two-dimensions  is 

P ■ P^expC-tinz  + i.kx  - -tut). 

For  a plane  wave  trav<  i.ng  down  the  z-axis,  this  becomes 

P ■ PQexp(^az  - -tut). 

The  transmitted  and  the  reflected  waves  together  form  a more  complete 
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solution  to  equation  (1); 

P ■ P+(x,*)exp(<tmz)  + P (x,z)exp(--£mz)  exp(-<j*»t). 

P+  designates  the  transmitted  wave  and  P the  reflected  wave.  Kote  that 
both  are  now  functions  of  location.  If  they  do  not  vary  rapidly,  more 
computer  space  may  be  saved  by  using  a coarse  grid  for  the  model  profile. 

The  transmitted  wave  is  separated  from  the  reflected  wave  and  nu- 
merical equations  are  developed  to  describe  its  behavior.  These  are 
solved  iteratively  with  a row-by-row  process  that  uses  the  left  and  right 
boundary  conditions  and  the  motion  found  on  the  previous  two  rows.  In 
this  way  the  wave  patterns  are  developed  over  the  entire  velocity  model. 
The  wave  amplitude  values  at  each  model  point  can  then  be  used  to  gener- 
ate the  reflected  waves  by  a similar  row-by-row  difference  procedure. 

Figure  6 and  Figure  7 were  plotted  using  the  procedure  for  trans- 
mitted waves.  Both  patterns  are  formed  by  a P-wave  moving  from  the  bottom 
to  the  top  of  the  illustrations.  The  scatterer  represents  a flat-floored 
chamber  filled  with  mud.  In  Figure  ft,  the  wavelengths  are  much  smaller 
than  the  dimensions  of  the  cavity.  The  wave  patterns  conform  quite  well 
with  ray  theory  and  the  focusing  of  this  mud  lens  is  clearly  shown.  In 
Figure  7,  the  radius  of  the  cave  is  smaller  than  the  length  of  the  wave 
and  diffractions  are  shown  to  be  an  important  feature.  A time  delay 
along  the  wave  front  is  also  present. 

Even  though  these  two  figures  Illustrate  refraction,  diffraction, 
and  interference  patterns,  there  is  a great  deal  of  information  that  has 
been  lost.  In  the  Claerbout  method,  shear  waves  are  neglected.  This 
prevents  the  formation  of  node  converted  waves  by  the  scatterer.  These 
S-waves  would  be  an  Important  factor  in  the  interference  patterns  found 
above  a cavity.  Also,  without  S-veves,  an  accurate  description  of  ground 
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Figure  6.  Wave  pattern  when  cavity's  radius  is  greater  than  wavelength. 
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Figure  7.  Wave  pattern  when  radius  of  cavity  is  less  than  wavelength. 
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motion  is  impossible.  Another  difficulty  with  the  Claerbout  method,  is 
that  due  to  the  approximations  mad*,  the  amplitudes  of  the  waves  become 
unstable  as  the  wave  deviates  from  vertical.  This  is  serious  since 
voids  often  develop  large-angle  diffracted  waves.  Even  with  the  large 
computer  costs,  the  best  way  to  study  scatter  is  with  a complete  finite 
difference  technique. 

As  a continuation  of  the  research  sponsored  by  the  Army,  the 
Ph.D.  dissertation  to  be  completed  by  Cavanaugh  will  be  on  the  use  of 
finite  difference  for  the  development  of  synthetic  seismograms  above 
air-,  water-,  and  mud-filled  cavities.  The  objectives  of  this  computer 
modelling  will  be:  (1)  to  obtain  a better  understanding  of  the  "turbu- 
lence" created  by  a wave  sweeping  past  a cavity;  (2)  to  find  if  it  is 
possible  to  delineate  the  shape,  size,  and  depth  of  a cavity  by  the  scat- 
ter that  it  causes;  (3)  to  determine  how  to  enhance  the  seismic  signature 
of  a cavity;  and  (4)  to  find  what  combination  of  size  and  depth  will  make 
the  disturbance  too  insignificant  to  be  recognized  on  a seismogram. 


CONCLUSIONS 


Biot's  work  demonstrated  the  existence  of  cylindrical  surface  waves, 
but  the  manner  of  their  propagation  makes  them  unusable  for  the  detection 
of  cavities.  Even  if  these  waves  were  recorded  above  a shallow  tunnel, 
their  velocity  and  particle  motion  would  make  them  indistinguishable 
from  the  Rayleigh  wave.  For  deep  tunnels,  they  would  not  be  noticed 
because  of  their  rapid  attenuation.  In  addition,  the  axial  surface  wave 
is  a travelling  wave  and  cannot  be  used  to  explain  resonance.  It  was 
thought  that  circumferential  surface  waves  on  a cylinder  could  explain 
these  oscillations,  but  they  cannot  be  shown  to  exist  theoretically . The 
long-lived,  periodic  vibrations  measured  above  some  cavities  appear  to  be 
due  to  flexural  waves.  These  can  only  be  formed  under  the  special  cond- 
tion  that  the  roof  of  the  cavity  behaves  as  a vibrating  plate,  thus  re- 
stricting their  occurrence  to  extremely  shallow  caves. 

For  caves  that  have  substantial  roofs,  the  anomalies  of  amplitudes 
and/or  arrival  times  that  are  occasionally  detected  in  a seismic  "shadow" 
zone,  may  be  a possible  means  of  delineation.  This  "shadow"  is  likely  a 
region  of  "turbulence"  created  by  the  cavity  acting  as  a scatterer  for  an 
up-coming  reflection.  Wave  scatter  cannot  be  adequately  investigated  by 
analytical  methods,  because  of  the  necessity  of  using  very  simple  boundary 
conditions.  The  inability  to  complete  the  problem  exactly  in  an  infinite 
half-space  makes  it  impossible  to  predict  the  interaction  of  the  scattered 
wave  with  a plane  surface.  Claerbout's  procedure,  which  is  an  approxima- 
tion to  the  finite  difference  method,  can  be  used  for  a model  containing 
many  Interfaces  with  complex  shapes,  but  its  solutions  can  become  unstable 
and  the  important  behavior  of  shear  waves  is  neglected.  To  obtain  the  best 


mathematical  description  of  cavity  wave  scatter,  it  will  be  necessary  to 
use  the  complete  finite  difference  method.  This  is  presently  being  done 
as  &n  extension  to  the  work  completed  with  the  support  of  the  Army. 

In  summary,  the  following  theoretical  approaches  have  been  explored 
in  this  study: 

1.  Resonance; 

a)  void  cavities  in  infinite  solids  (intrinsic  cavity  resonance), 

b)  void  cavities  in  infinite  half-spaces  (resonance  caused  by  the 

interaction  between  cavity  and  land  surface), 

c)  axial  cavity  surface  waves  observable  at  land  surface, 

d)  circumferential  cavity  surface  waves  observable  at  the  land 
surface, 

e)  the  cavity  roof  acting  as  a vibrating  mebrane,  and 

f)  the  cavity  roof  acting  as  a vibrating  thin  plate. 

2.  Scattering; 

a)  analysis  by  ray  tracing  using  Snell's  Law  in  three-dimensions  and 

b)  analysis  by  finite  difference  approximation  to  the  wave  equation. 
Only  the  thin  plate  model  and  the  finite  difference  approach  seem  to  hold 
promise.  The  thin  plate  model  would  apply  only  to  cavities  very  near  the 
surface.  For  deeper  caves,  the  finite  difference  approximation  for  wave 
scatter  offers  the  best  possibility  of  producing  a general  method  of  de- 
tection. This  has  yet  to  be  fully  developed. 


INTRODUCTION 


The  experimental  research  conducted  under  this  ('.rant  was  directed 
towards  the  development  of  a seismic  field  procedure  for  the  detection 
and  delineation  of  subterranean  cavities.  The  basic  underlying 
presumption  of  this  study  was  that  a resonance  phenomenon  is  associated 
with  voids  when  they  are  subjected  to  an  induced  seismic  field.  This 
presumption  was  based  on  experimental  evidence  reported  by  Watkins  et  al 
(1967),  and  by  Godson  and  Watkins  (1968),  and  on  theoretical  evidence 
presented  by  Stewart  (1971). 

To  illustrate  the  resonance  phenomenon  associated  with  voids,  a 
seismogram  over  a lava  tunnel,  as  reported  by  Watkins  et  al^  is  reproduced 
in  Figure  1,  They  reported  that  the  high  amplitude  oscillations,  as 
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Figure  I.  from  Watkins  et  al  (1967) 
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shown,  were  recorded  most  intensely  directly  over  the  cavity,  and  that 
this  result  was  similar  to  the  seismograms  recorded  over  two  other 
near-surface  cavities  that  were  studied.  They  noted  that  the  dominant 
frequency  is  relative  uniform  with  tine,  and  that  the  in-phascncss  of 
the  oscillations  was  indicative  of  resonance. 

Figure  1 in  itself  presents  a very  strong  argument  for  cavity 
resonance,  and  based  on  this  result  the  initial  research  effort  was 
directed  towards  a more  accurate  definition  of  this  phenomenon. 
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EXPERIMENTAL  PROCEDURE 


Test  Sites  The  bulk  of  the  seismic  investigations  were  conducted 
over  segments  of  Cathedral  Caverns  and  Onondaga  Cave,  located  approximately 
seven  miles  south  of  Leasburg,  Missouri.  Measurements  were  also  taken 
at  three  other  test  sites  over  known  cavities,  cave  systems,  or  suspected 
cave  systems,  but  have  yielded  little  results.  In  these  areas  the 
principal  problem  was  the  uncertainity  of  the  siesmic  traverse  relative 
to  the  underground  void,  excessive  depth  of  the  cavity  system,  and  the 
cavity  size.  Consequently,  results  will  only  be  presented  for  the 
Cathedral -Onondaga  areas. 

A plan  view  of  the  Cathedral -Onondaga  systems  is  shown  in  Figure  2. 

Seismic  lines  that  shall  be  discussed  in  this  report  are  shown  in  the 
figure  as  dark  lines  designated  as  Tl,  T2,  and  T3. 

Traverse  T2  is  over  a known,  and  accurately  surveyed,  cave  passage. 

The  trace  of  this  passage  above  ground  was  determined  by  means  of  surveying. 

The  passage  in  this  segment  of  the  tunnel  is  rather  unifoim  in  shape  with 
a height  of  roughly  10  feet  and  a width  of  51  feet.  The  roof  is 
located  approximately  162  feet  below  the  surface. 

Traverse  Tl  was  seismically  surveyed  to  detect  a possible  extension 
of  the  western  segment  of  Cathedral  Caverns.  Traverse  T3  was 
conducted  in  the  vicinity  of  Oi'iondaga  Cave  in  an  attempt  to  detect  an 
extension  of  a water  filled  segment  of  the  cave.  Both  Tl  and  T3  were 
traverses  over  previously  unexplored  areas. 

During  the  course  of  the  three  years  of  this  grant,  and  for  two  years 
prior  to  the  initiation  of  this  grant,  these  test  sites  were  repeatedly 
investigated  many  times,  with  slight  modifications  of  the  traverse  orientation, 
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improved  instrumentation,  or  principally  to  test  a reorientation  of 
the  research  direction. 

Instrumentation  Ground  velocity  was  measured  with  Sprengnether, 

Model  S-6000,  three -component , short  period  seismometers,  and  digital 
grade  vertical  and  horizontal  seismometers.  The  seismometer  outputs 
were  through  programmable  gain  amplifiers,  processed  through  a 

32  channel  analog  to  digital  converter  and  NOVA  800  Jumbo  Computer, 
and  permanently  stored  on  a moving  head  Diablo  cartridge  disk  system. 

The  data  acquisition  system  and  reduction  system  is  described 
in  Appendix  A of  this  report.  The  acqusition  system  is  housed  in  a 
seismic  van  as  shown  in  Figure  3.  The  basic  computer  acquisition  system 
as  installed  in  the  van  is  shown  in  Figure  4. 

The  frequency  response  of  the  seismic  system  ?s  Hat  from  2 to  200 
Hertz  for  the  system  using  Sprengnether  seismometers.  Using  the 
digital  grade  seismometers  the  response  of  the  system  is  flat  over  the 
range  of  8 to  200  Hertz. 

In  the  results  presented  a chemical  explosive  was  used  as  the  energy 
source.  Early  investigations  employed  a mechanical  vibrator. 

However  the  energy  production  of  the  instrument  was  far  too  low  to  be 
useful . 


Fit’ure  . 


Seismic  van 


RESULTS 


The  cavity  underlying  Traverse  T2  was  chosen  as  the  principal 
scattering  object  to  be  studied.  This  cavity  exhibited  mud-free 
walls,  uniformity  in  shape  and  direction  over  a large  portion  of  its 
length,  and  accurate  spatial  measurement  control.  Moreover,  this 
roughly  cl iptic-cyl indrical  void  was  not  too  large,  nor  too  small,  and 
was  located  at  an  interesting  depth  of  162  feet  below  the  surface. 

Cavities  of  this  size  and  depth  have  hitherto  been  opaque  to  geophysical 
probes . 

The  principal  results  of  this  investigation  were  obtained  over  this 
tunnel,  after  many,  many  trials.  These  results  are  shown  in  Figure 
5 which  presents  portions  of  the  three -component  measurements  along 
Traverse  T2.  The  measurements  shown  were  taken  at  10  foot  intervals 
along  a line  oriented  at  a right  angle  to  the  cavity  axis.  Seismometer 
#6  was  located  directly  over  the  center  of  the  cavity.  The  shot 
point  was  50  feet  from  Station  1. 

From  the  volume  of  data  taken  in  the  near  proximity  of  Traverse  T2, 
many  different  seismic  indications  of  the  presence  of  a cavity  were 
discovered.  Some  of  these  indicators  are  discussed  in  Rechtien  and 

Stewart  (1975),  which  include  diroct  scattered  P and  S wave  energy  from 
the  cavity  roof,  and  the  predominance  of  horizontal  notion  in  planes 
parallel  to  the  tangent  planes  to  the  cavity  at  late  arrival  times. 

These  results  still  stand.  However  the  analysis  of  early  arrival  times 
is  subject  to  considerable  error,  particularly  in  this  area  where  the 
topography  is  rugged  and  weathering  depths  and  velocity  vary  considerable 
over  short  distances. 
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On  the  other  hand,  the  horizontal  motion  in  planes  parallol  to 
the  cavity  boundary  have  been  consistently  observed,  often  witii  energies 
far  exceeding  the  SH  wave  energy  proceeding  from  the  sourco.  In 
addition  a reflected  Love  wave  from  cavities  has  been  observed,  but 
not  consistently.  The  only  type  of  seismic  o/ent  that  has  been 
measured  consistently  is  illustrated  in  Figure  5.  hut  before 
describing  these  results  it  might  be  worth  stating  that  ovor  all  air 
filled  cavities  with  mud  free  walls  that  were  investigated,  never  was 
a record  obtained  that  exhibited  the  "resonance"  signature  of  Figure  1, 
as  reported  by  Watkins  et_  al  . However  the  cavities  in  the  area  of 
T2,  and  at  other  locations,  were  considerably  doeper  than  those  reported 
be  Watkins  e£  al_.  And  perhaps  the  shallowness  of  the  cavity  that 
they  investigated  is  a significant  factor  to  their  results. 

Figure  5 presents  the  results  of  three-component  ground  velocities 
from  400  to  700  milli-seconds.  Earlier  data  is  predominated  by 
surface  wave  energy  which  shows  normal  transmission  behavior.  Later 
data  shows  random,  scattered  plurse  arrivals  and  is  of  little  interest. 

The  data  as  shown  has  boon  processed  through  a digital  band  pass 
filter  after  the  data  was  windowed  by  a function  which  allotted  suppression 
of  the  data  below  350  milli-seconds  and  above  l.'lOO  milli-seconds.  The 
windowing  was  performed  to  suppress  the  effect  of  the  high  energy  surface) 
waves  on  the  band  pass  filtering  process.  The  inclusion  of  the  high 
frequency  data,  above  35  Hertz,  seriously  distorts  the  data.  Although 
the  basic  pattern  illustrated  in  Figure  5 is  there,  it  is  not  at  all 
obvious  to  the  naked  eye. 
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The  first  indication  of  a cavity  is  shown  on  the  Vertical  data 
i’.v  events  that  arc  concave  upward  in  the  central  portion  (500-000  ms) 
of  the  seismogram.  Note  the  negative  (un- shaded)  event  which  begins 
at  Station  1 at  550+  ms  and  ends  at  Station  11  at  the  same  tijnc  and 
dips  to  580+  ms  at  Station  o.  The  following  negative  event  parallels 
this  some  85  ms  later,  and  the  negative  event  preceding  shows  the 
same  concavity  upwards  on  the  rlanu..  o*'  the  cavity,  ns;t  in  fact  it  in 
convex  downward  over  the  central  portion  of  the  void  at  Stations  6 and  7. 
A dead  spot  results  beneath  this  convexity  where  the  ]>ositivc  (shaded) 
energy  diminishes  to  near  zero  through  destructive  interference.  The 
downward  convexity  of  the  positive  event  directly  above  this  dead  spot 
is  even  more  pronounced.  Also  note  t.n  nigh  energy  events  directly 
above  this  shadow  zone  between  Stations  5 through  3 between  425  to  47f>  ms 
This  high  energy  event  shows  no  phase  relationship  with  the  traces  on 
either  side. 

On  the  longitudinal  data  it  is  obvious  that  Stations  5 and  8 are 
points  of  phase  reversals  across  the  records  between  450  to  600  ns. 

At  640  ms  Stations  1 through  4 and  0 through  11  have  positive  (shaded) 
phase,  whereas  Stations  6 and  7 arc  reversed  in  phase,  and  Stations 
5 and  8 are  somewhat  nulled.  This  behavior  can  also  be  followed 
at  earlier  times  to  450  ms. 

The  Transverse  data  show  little  or  nothing  in  this  record  interval. 


If  anything  it  shows  a diminishing  of  wave  energy  in  the  central  portion 
of  the  traverse. 


Particle  velocity  hodograms  for  Stations  3,  6,  and  9 for  the 
seismic  interval  of  520  to  600  milliseconds  are  shown  in  Figure  6. 

Those  hodograms  are  typical  of  the  entire  traverse  and  seem  to 
imply  that  t lie  vertical  and  longitudinal  motions  are  independent. 

This  conclusion  is  not  too  surprising  since  the  Wave  liquation  in 
a cylindrical  coordinate  system  naturally  separates  motion  in  the 
r-0  plane  to  he  independent  of  axial  motion.  Consequently  the 
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results  are  more  meaningful  as  viewed . from  the: time  traces  themselves. 

The  identification  of  a diagnostic  event,  or  sequence  of  events, 
as  given  in  Figure  5,  is  the  principal  result  of  this  study. 

Although  the  mechanism  that  results  in  this  cavity  signature  on  a 
seismogram  has  not  been  identified,  it  seems  probable  that  it  is  a 
surface  wave  phenomenon,  either  associated  with  cavity  interference 
with  a deep  penetrating  surface  wave  travelling  away  from  the  source, 
or  from  surface  waves  associated  with  the  cavity  itself  and  generated 
by  mode  conversion  from  waves  reflected  from  doep  horizons,  or  from 
entrapment  of  passing  surface  wave  energy. 

The  data  of  Figure  7 siicds  considerable  light  on  this  question, 
and  very  clearly  supports  the  deep  penetrating  surface  wave  interference 
line  of  thinking.  This  seismogram  was  obtained  over  a segment  of 
Cathedral  Caverns  in  the  vicinity  of  Travorsc  Tl.  Fifteen  seismic 
stations  are  shown  and  one  three-component  measurement  which  was  taken 
at  a point  near  Station  1.  This  record  has  not  been  filtered,  and  in 
fact  the  frequency  spectra  of  t.u-s,-  r.vacos  are  siown  in  Fi^uv..  . 

This  particular  shot  was  not  in-line  with  the  traverse. 


Figure  1.  SeisBOgre*  of  vertical  Motion,  Shot  1,  line  U. 
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Figure  8.  Frequency  Spectra,  Shot  1,  Line  11. 


? 


61 


In  reference  to  Figure  7,  a dead  S|X5t  is  shown  in  the  vicinity  of 
Station  6 around  600  ms.  Station  ft  overlies  a cavity  of  approximately 
the  same  site  and  depth  of  burial  as  that  of  Traverse  T2.  It  seems 
presumptuous  to  idontify  a wave  on  this  record  by  drawing  a straight 
line,  but  for  the  purpose  of  focusing  on  a particular  segment  of  the 
record  a straight  line  has  been  drawn.  This  straight  line  goes 
through  points  of  common  phase  on  lioth  sides  of  Station  C,  although 
the  data  to  one  side  of  Station  6 is  18')  degrees  out  of  phase  with 
that  on  the  other  side.  If  one  would  calculate  velocities  from  a 
delta-t  - delta-x  analysis,  an  extremely  slow  velocity,  in  the  range 
of  7i:0  feet  per  second  would  be  obtained.  Yet,  the  three-dimensional 
counterpart  of  the  680  ms  wiggle  on  trace  one  clearly  exhibits  the 
90  degree  vertical -horizontal  phase  shift  of  a Rayleigh  wave. 

Moreover,  the  data  is  clearly  indicative  of  waves  travelling  to  the 
right,  away  from  the  source,  as  they  should.  It  appears,  therefore, 
that  the  dead  spot,  at  least  in  this  instance,  is  due  to  the  interference 
cf  the  cavity  with  a passing,  deep  penetrating  surface  wave.  It  is 
hard,  however,  to  account  for  such  a slow  speed.  Moreover,  the 
concave  and  convex  curvatures  and  high  energy  events  above  the  void, 
which  were  clearly  evident  on  the  data  of  Traverse  2,  are  not  present  in 
this  record. 

In  reference  to  Figure  8,  the  surface  wave  energy  centered  at 
approximately  16  Hertz  is  clearly  diminished  over  the  cavity  This 
result  was  common  to  the  data  of  Traverse  2 as  well  as  to  all  the  data 


taken  over  knc>n  cavities. 
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LINE  A 

Figure  y.  Gravity  profile  over  Traverse  T2 
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Kefore  leaving  Traverse  2 it  in  instruct? ve  to  consider  the 
gravity  anonaly  over  this  traverse.  This  data  Is  shown  in  Figure 
9,  with  Gravity  Station  21  corresponding  to  Seismic  Station  1,  and 

Gravity  Station  11  corresponding  to  Seismic  Station  11.  There 

is  0,  faint  gravity  anomaly  in  the  vicinity  of  Gravity  Station  13. 

This  nnoawly  is  hardly  not i cub ^ in  tiie  observed  or  free-air  gravity 

profile,  and  in  fact  one  wonders  if  the  anomaly  is  not  just  a correction 

effect.  Regardless,  grsvity  hardly  s«ems  like  a feasible  tool 
for  cavity  detection  of  this  size  and  depth.  These  measurements 
were  tfken  with  a Worden  Gravity  Motor. 

The  elevation  of  the  traverse  illustrated  in  this  figure  shows 
the  typical  problem  encountered  iu  the  search  for  cavities;  they 
are  usually  found  beneefh  steep  slopes  of  hillsides.  This  fact 
creates  serious  problems. in  the  interpretation  of  the  results. 

Figure  lb  stows  the  seismic  results  along  Traverse  T1 . This 
data  was  actually  acquired  prior  to  the  initiation  of  this  grant  over 
an  area  of  an  unexplored  region  of  Cathedral  Caverns.  Itocaivse  of 
the  numerous  indications  of  cavities  along  this  traverse,  this  area 
was  chosen  to  test  t series  of  broadside  shots  in  ;»n  effort  to  obtain 
three-dimensional  coverage  of  th?  subsurface  using  CUP  techniques. 

In  Figure  1G,  surface  waves  dominate  the  record  down  to  ab>ut 
700  cs  in  both  the  forward  and  rnverse  shots.  On  the  North  shot 
data  a dead  spot  first  occurs  near  Station  7 at  around  250  ms.  and 
again  in  the  vicinity  of  Station  12  ut  500  ms,  and  at  Station  25  at 
700  ms.  These  all  appear  to  bo  disruptions  of  the  passing  surface 

waves.  Station  25  on  tlio  South  shot  also  nppeers  to  be  in  the 
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neighborhood  of  a dead  spot  at  about  600  ms,  and  Station  21  at  250  nis 
and  1050  ms.  One  can  also  see  many  areas  of  convex  and  concave 
curvature  in  the  records.  Two  clear  examples  can  be  seen  on  the 
South  shot;  the  convex  curvature  between  stations  1‘)  and  27  at  about 
570  ms,  and  the  concave  curvature  underlying  the  dead  spot  at  700  ms; 
then  again  the  convexity  at  1000  ms  between  Stations  17  and  24,  and 
the  underlying  concavity  at  1100  ms.  It  is  very  clear  that  at  least 
the  tops  of  these  dead  spots  are  the  terminal  end  of  a passing  surface 
wave.  These  data  are  much  like  the  data  of  Traverse  2,  except  for 
the  high  energy,  inphase  events  directly  over  the  cavity.  Recent 
exploration  of  Cathedral  Caverns  has  resulted  in  a subterranean  passage 
in  this  area,  although  its  precise  location  is  not  known. 

Very  strong  phase  reversals  on  traces  7 and  14  in  the  neighborhood 
of  1000  ms  on  the  North  shot  data. can  be  seen.  This  reversal  can 
also  be  observed  at  Station  8 at  800  ms  on  the  South  shot  data.  This 
behavior  is  similar  to  the  results  of  Figure  7.  The  rather  strong 
reverberations  in  this  region  of  the  North  record  and  the  reverberations 
between  stations  7 and  11  on  the  South  shot  is  somewhat  indicative  of 
the  resonance  reported  by  Watkins  ct  al_.  There  are  no  known  cavities 
in  this  region.  liowever,  a gravity  survey  over  this  traverse  was 
performed,  and  the  results  are  shown  in  Figure  11.  The  gravity 
and  seismic  stations  in  this  case  coincide,  and  the  data  show  a gravity 
low  both  otar  the  area  of  reverberation  near  Station  9 and  near  Station  21. 
The  question,  "Why  the  reverberation?"  arises.  A cavity  signature 
meets  ail  the  criteria  near  st.it  ion  21.  but  what  is  happenin',  near 
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Figure  il.  Gravity  profile  ov«r  Traverse  Tl. 
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The  data  from  Traverse'  T3  answers  this  question.  This  data 
is  presented  in  Figure  12.  In  reference  to  this  figure  a duplicate 
of  the  results  of  Watkins  et  al  is  seen.  This  profile  was  taken 
over  an  area  suspected  of  cavities,  and  upon  the  attainment  of  this 
result  a hole  was  drilled  at  Station  13.  The  drill  penetrated  15 
feet  of  rock  before  sinking  120  feet  into  cave  mud  before  drilling 
was  terminated.  A solid  bottom  was  not  encountered. 

Gravity  and  electrical  resistivity  profiles  were  obtained  over 
this  profile  and  are  presented  in  Figure  13.  The  position  of  the 
drill  hole  is  indicated  at  Gravity  Station  34.  This  area  is  both 
a gravity  and  electrical  resistivity  low,  well  indicative  of  mud. 

Parallel  resistivity  profiles  were  conducted  in  an  effort  to  determine 
if  a fault  plane  or  a vertical  filled  solution  cavity  had  been 
encountered.  The  pattern  of  measurements  is  shown  in  Figure  * , 
and  the  results  are  shown  in  Figures  15  and  16.  These  resul  show 
the  circular  nature  of  the  mud  fill,  which  points  to  a cylindrical 
solution  cavity. 

From  the  results  of  Traverse  T3  one  questions  if  the  cavity  studied 
by  Watkins  et_  al_  from  which  they  obtained  the  data  of  Figure  1 was  indeed 
empty  as  they  claim.  The  data  from  T3  is  almost  identical  to  theirs. 

From  this  result  it  appears  resonable  to  assumo  that  if  strong  reverberations 
are  recorded,  such  as  those  exhibited  by  Figure  10  for  Traverse  T2,  that 
a mud  or  water  filled  void  has  been  encountered.  This  is  the  only 
case  for  which  "resonance"  is  physically  plausible  anyway. 
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Figure  12.  Seismic  data  over  Traverse  T3.  Spacing 
10  foot  intervals.  Timing  interval  100  ms.  Shot 
point  50  feet  from  end  of  line  at  Station  1. 
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PLAN  VIEW  OF  GRAVITY  AND  RESISTIVITY  TRAVERSE 
INTERSECTIONS  - LOST  CAVE  RIDGE 

Figure  14.  Plan  view  of  Electrical  Resistivity  sapping. 
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ELECTRICAL  RESISTIVITY 
LIMES  A,  B,  C,  D,  E,  F,  C,  - LOST  CAVE  RIDGE 
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Figure  IS.  North-south  Electrical  Resistivity  profiles. 
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Figure  16.  Bast -west  Electrical  Resistivity  profiles. 
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CONCLUSIONS 

The  data  presented  show  two  very  clear  cut  examples  of  cavity 
signatures  in  seismic  data;  one  associated  with  a void  cavity,  and 
one  for  a cavity  filled  wir.h  mud.  These  signatures,  moreover, 
nave  been  observed  on  other  traverses,  over  known  or  suspected 
cavity  areas  where  the  presence  of  cavities  have  subsequently 
been  verified  by  subsurface  tunnel  exploration  or  through  gravity 
measurements. 

Tne  seismic  manifestation  of  a cavity  on  a seismogram  appears 
to  lie  the  re- ult  cf  the  interference'  of  the  cavity  with  passing 
surface  waves.  It  can  be  observed  in  the  data  that  the  early 
Rayleigh  and  love  waves  associated  with  shallow  layers  that  overlie 
the  cavity  propagate  relatively  undisturbed.  The  deeper 
penetrating,  later  arriving  surface  wave  energy,  on  the  other  hand, 
very  clearly  exhibits  the  interference  effect  of  the  void.  The 
very  low  frequency  of  the  seismic  events  manifesting  the  cavity  also 
point  to  surface  wave  energy.  Moreover,  the  abruptness  of  the 
lateral  variations  of  the  data  imply  that  the  seismogram  is  reflecting 
the  very  near  surface  velocity  and  density  variations,  as  opposed  to 
the  deeper  structural  variations. 

It  was  initially  believed  that  the  cavity  signature  was  a result 
of  the  interference  of  the  void  with  a reflected  wave  returning  to 
the  surface  fro«  the  basement  complex,  which  for  this  area  is  less 
than  3,000  feet  in  depth.  With  an  assumed  average  section  velocity 
of  as  low  as  6,000  feet  per  second  the  reflected  events  would  have  to 
occur  at  record  times  of  less  than  one  second.  However,  the 
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velocity  is  probably  higher  than  6,000  feet  per  second  and  it*s  hard 
to  account  for  the  later  events  on  the  record  in  terms  of  reflected 
wave  interference.  Moreover,  from  such  a deep  reflector,  relative 
to  the  spread  length  of  the  seismic  line  used  in  this  study,  the 
north  and  south  shots  should  show  a similarity  of  the  reflected 
events  as  they  are  distorted  in  passing  shallow  vvloeutj' 

i rreo.; I ari  t ies . The  ray  paths  wo  ilii  ever,  J>c  near  vortical  . 

The  dissimilarity  of  the  records,  the  low  frequency  of  the 
cavity  signature,  and  the  very  abrupt  lateral  variations  of  the 
events  all  but  negate  the  plausibility  of  a reflected  wave  phenomenon. 
The  surface  wave  interference  theory  seems  much  more  reasonable. 

The  data  show  that  cavities  as  small  as  10  feet  in  height, 

SI  feet  in  width,  and  162  feet  below  the  surface  can  be  detected 
seismically.  However,  the  corresponding  seismic  signature  is 
not  all  that  obvious  in  a reeord  section.  For  example,  a very 
complete  set  of  seismic  data  obtained  with  broadside  shooting 
using  common  depth  point  coverage,  12  fold,  was  acquired  along 
a traverse  100  feet  froa,  and  paralleling, Traverse  T1 . But  the 
data  was  not  presontad  in  this  report  because  the  data  did  not 
show  any  obvious  indication  of  the  cavity  which  is  known  to  exist 
beneath  this  traverse.  Only  one  record,  which  is  presented  in 
Figure  7 of  this  report,  showed  the  presence  of  the  cavity,  but  it 
also  didnot  exhibit  the  typical  cavity  signature.  This  reeord 
was  from  a shot  furthest  removed  from  the  iinc,  anJ  represented 
the  most  in-line  configuration. 

This  result  suggests  that  two  criteria  are  important  in  acquiring 
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seismic  data  over  cavities.  First,  the  shot  point  and  receivers 
must  be  in-line.  Secondly,  perhaps  the  seismic  line  must  lie 
at  a right  angle  to  tho  horizontal  cavity  axis.  The  first 
criterion  can  readily  be  met  in  the  field.  The  second  is  impossible 
to  achieve  in  unexplored  area:-.. 

Cavity  detection  by  seismic  means  is  possible.  Yet  the  procedure 
is  tedious.  The  required  close  spacing  of  the  seismometer 
stations  prohibits  its  use  as  a reconnaissance  tool.  The  method 
could  be  applied  effectively  over  small  areas,  for  instance  a dam 
site,  lagoon  site,  or  the  like.  Out  the  volume  of  data  required 


for  large  areas  would  be  prohibitive. 
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MINICOMPUTERS  IN  THE  ACQUISITION  AND  ANALYSIS  OP  SEISMIC  DATA 

J.  H.  Tracey,  R.  D.  Rachtian  and  K.  L.  Hambacker 
Univaraity  of  Missouri -Roll a 
Holla,  Missouri  (5401 


INTRODUCTION 

This  paper  describes  the  use  of  a mobile 
laboratory  and  minicomputer  ays terns  to  sup- 
port geophysics  research  and  instruction  at 
the  University  of  Missouri-Rolla.  The  equip- 
ment includes  a mobile  van,  a field  data  ac- 
quisition system  and  a minicomputer  network 
for  data  analysis  and  classroom  replay. 

SEISMIC  STUDIES 
Generalised  Seismic  Process 

Seismic  methods  sre  employed  to  probe 
the  internal  stratification  of  the  earth  in 
an  effort,  generally,  to  detect  and  isolate 
abnormal,  localised,  geological  structures 
that  may  exist  near  the  surface1.  Such 
structures  might  be  fault  zones,  solution 
cavities,  salt  domes,  buried  stream  channels, 
or  perhaps  simply  undulations  in  the  rock 
layers  which  might  be  favorable  to  the  ac- 
cumulation of  oil  or  mineral  deposits.  Engi- 
neering problems,  such  as  the  response  of 
tall  buildings  or  dams,  to  seismic  energy, 
also  fall  within  the  broader  scope  of  seis- 
mology. Thus,  seismic  methods  deal  with  the 
transmission  of  elastic  waves  through  the 
heterogeneous  earth  and,  occasionally,  their 
effects  on  man-made  structures.  * 

For  the  application  of  seismic  methods 
to  non-conventionsl  subsurface  problems,  it 
is  essential  that  the  experimental  test  pro- 
cedures be  continually  modified  upon  the  ac- 
quisition of  new  information,  in  order  to  ac- 
complish a given  research  task.  The  research- 
er must  have  the  flexibility  to  sero-in  on 
the  localised  area  of  interest,  and  thus  avoid 
the  tremendous  volume  of  data  required  for 
blind  wider -rat, gs  ground  coverage.  This 
flexibility  requires  that  a decision  making 
capability  be  made  an  integral  pert  of  the 
test,  and  as  such  requires  an  immediate  data 
analysis  capability  upon  acquisition.  The 
needed  system  cepabilities  for  near-surface 
exploration  are<  the  hodograph,  or  particle 
trajectory  plots,  used  prim’arily  for  wave 
type  identification)  the  Fast  Fourier  Trans- 


form, used  principally  for  spectral  analysis, 
band-pass  filtering,  notch  filtering,  or 
spectral  shaping)  cross-correlation)  decon- 
volution) and  optimum  filtering,  used  for 
shaping  the  frequency  spectrum  according  to 
some  knowledge  of  prior  filtering  actions  or 
according  to  some  assumptions  about  signal 
characteristics  or  desired  output.  These 
features  provide  a wide  range  of  field  analy- 
sis capability. 

DESCRIPTION  OF  MOBILE  LABORATORY 

The  mobile  laboratory  consists  of  a con- . 
verted  2 1/2  ton  U.S.  Army  truck  with  a com- 
pletely enclosed  rear  compartment.  The  van' 
houses  the  minicomputer,  its  peripherals,  and 
associated  test  equipment.  The  van  is  com- 
pletely air  conditioned  and  capable  of  being 
powered  from  a 220-volt  utility  outlet  or  from 
an  on-board  generator. 

The  minicomputer  is  a NOVA  800  Jumbo  Com- 
puter with  16K  core,  and  hardware  multiply  and 
divide  option.  A moving  head  Diablo  cartridge 
disk  system  with  1.2  million,  16-bit  word 
capacity  and  a Tektronix  4010  storage-tube 
graphics  terminal  are  peripheral  to  the  NOVA. 

A high-speed  reader  is  included  for  initial 
loading  of  application  software  and  utility 
routines.  The  system  includes  an  analog-to- 
digital  converter  which  is  capable  of  convert- 
ing up  to  eight  analog  channels  at  a maximum 
rate  of  30  KHz  for  a single  channel. 

A Texas  Instruments  DSP  (Digital  Field 
System)  data  recorder  is  also  available  in  the 
van.  This  unit  permits  the  digitization  of  up 
to  32  analog  channels  which  can  later  be  pley- 
ed  back  into  the  NOVA.  The  DFS  unit  is  used 
as  a back-up  system  for  the  NOVA  and  is  also 
ujed  in  those  experiments  requiring  the  con- 
tinuous digitisation  of  more  than  eight  chan- 
nels. 

FIRED  DATA  ACQUISITION  AMD  ANALYSIS  SYSTEM 
Software  Development 

The  software  written  for  the  min. computer3 
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was  designed  to  Mat  a number  of  system  con- 
•traints  in  terms  of  program  alia,  a pa ad  of 
execution,  aaaa  of  oparation  and  invulner- 
ability  to  uaar  arror  aa  wall  aa  aafaty  of 
storad  data.  Tha  program  ia  complataly  core 
raaidant  along  with  a copy  of  tha  Data 
Ganaral  Corporation  Mai  Tima  Oparating  Sys- 
tam  for  a total  langth  of  undar  IX  word a of 
memory.  Thia  laavaa  IK  words  of  cora  for 
data  a tor  age,  a tranaform  array,  dials  I/O 
buffer,  and  command  sequence  atoraga.  Spaed 
of  execution  received  high  priority  in  tha 
software  design  to  minimi is  uaar  annoyance. 
This  goal  of  instantaneous  machine  response 
is  a major  factor  In  keeping  tha  entire  pro- 
gram core  resident.  Tha  command  structure 
is  designed  to  minimise  keyboard  entries  for 
maximum  user  time  efficiency  at  the  cost  of 
a longer  user  familiarisation  period.  Single 
letter  mnemonics  are  used  to  initially  branoh 
to  the  various  areas  of  program  control  with 
all  further  information  being  prompted.  Con- 
siderable effort  has  been  expanded  to  make 
the  system  and  data  resistant  to  user  errors 
by  keying  sensitive  system  parameters  which 
control  buffer  placement  and  disk  write 
pointers.  As  long  as  this  software  key  is 
not  used  the  user  may  randomly  type  at  the 
keyboard  without  damaging  the  program  or 
stored  disk  data.  A break  key  is  defined 
which  allows  any  operation  to  be  interrupted 
with  a single  keystroke.  Since  considerable 
tiM  and  cost  is  involved  in  obtaining  data, 
reliable  storage  is  imperative.  For  this 
reason  data  is  stored  with  distributed  iden- 
tification information  so  that  damage  to  a 
critical  index  seotor  cannot  causa  the  loss 
of  a major  portion  of  the  disk  information. 
The  limited  amount  of  code  devoted  to  disk 
I/O  also  lowers  the  chances  of  a malfunction- 
ing program  or  maohine  getting  stuck  in  a 
loop  which  damages  storad  data  on  disk. 

System  Operation 

Bursts  of  data  from  up  to  eight  channels 
are  obtained  from  the  A/D  converter  at  ac- 
curately timed  rates  of  up  to  1 XHs  and  stor- 
ed in  oore  memory.  After  sampling,  a oopy  of 
the  raw  data  is  immediately  stored  on  disk 
while  the  data  in  oore  may  be  manipulated  for 
analysis,  individual  channels  of  data  may 


have  a DC  level  removed,  be  shifted  in  time, 
and  be  normalised  in  amplitude.  Channels  may 
be  added,  subtracted,  and  multiplied  col- 
lectively for  such  operations  as  signal  aver- 
aging. Arrays  of  time  samples  may  ha  normal- 
ised and  Fast  Fourier  Transformed  at  2*® 
samples  in  about  one  and  one  half  seconds. 

The  spectrum  may  be  modified  and  an  inverse 
transform  obtained.  Also,  single  events  of 
interest  in  a waveform  may  be  isolated  or 
used  as  models  for  wavelets  entered  from  a 
digitising  tablet. 

A variety  of  displays  are  available  of 
the  data  or  functions  of  the  data  without  any 
modification  to  the  data  arrays  used  or 
storage  of  the  resultant  functions.  All  plots 
may  be  scaled  by  powers  of  two  and  offset  to 
display  any  desired  window  of  information. 
Vertioal  and  horisontal  markers  are  scaled  to 
be  consistent  with  the  displays  for  waveform 
measurement  .*  Thous ind-point  clots  typically 
take  slightly  over  two  seconds.  Single  chan- 
nel time  plots  are  available  of  information 
stored  in  the  sample  buffer  and  of  real  infor- 
mation stored  in  the  transform  array.  If  the 
information  in  the  transform  array  is  rectan- 
gular complex  spectral  information,  magnitude 
and  phase  angle  functions  may  he  plotted  as 
well  as  the  stored  real  and  imaginary  repre- 
sentations of  the  transform.  Several  plots 
involving  more  than  one  channel  of  sampled 
Information  are  available,  windowed  plots  of 
particle  position  with  time  may  be  derived 
from  three  channels  of  data.  Oblique  and  pro- 
jected views  from  any  angle  are  allowed  at 
selectable  speeds  up  to  maximum  terminal 
plotting  speed.  Particle  plots  may  be  inter- 
rupted for  point  number  identification  and 
then  continued  to  allow  isolation  of  individ- 
ual events  in  the  three  dimensional  repre- 
sentation. Displays  are  also  available  of 
correlation  and  convolution  functions  of  any 
two  channels  of  sample  data. 

The  command  structure  of  the  program  al- 
lows rapid  and  convenient  use  of  various 
analysis  methods  and  hence  allows  the  user  to 
concentrate  on  innovative  techniques  for  data 
reduction  which  might  be  neglected  or  missed 
in  a system  with  slower  turnaround.  Once  a 
promising  sequence  of  analysis  steps  has  been 
arrived  at  it  may  be  stored  and  recalled  at 
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will  for  us*  on  other  data  aata.  This  faatura 
of  tha  program  has  provad  indiapanaabla  for 
tha  ordarly  display  and  analysia  of  three-axis 
data  where  tha  asms  oparationa  ara  required 
for  aach  channai.  Tha  command  aaquancaa  ara 
not  intended  for  uaa  in  building  baaio  new 
functions,  which  ara  baat  implamantad  in  as- 
samblar  coda,  but  maans  ara  aval labia  to  antar 
run  tima  variables,  do  aimpla  arithmatic,  do 
an  arithmetic  branch,  pause  for  displays  and 
sat  up  loops.  Thssa  command  aaquancaa  are 
stored  and  executed,  as  they  ara  initially 
entered,  a?  errors  ara  detected  immediately 
and  preparation  time  is  kept  to  that  normally 
raquirad  for  command  keyboard  entry. 

ROUS  Of  OMK  COMPUTER  NETWORK 

Figure  1 shows  a portion  of  tha  currant 
configuration  of  tha  UHR  minicomputer  net- 
work3 which  complements  the  data  acquisition 
system  described  in  tha  previous  saotions. 

Tha  network  consists  of  a NOVA  minicomputer 
co-locatad  with  and  attached  to  tha  multi- 
plexor channai  of  an  IBK  360/S0  Computer. 
Several  remote  minicomputers  ara  attached  via 
high-speed  lines  (50K  baud)  to  this  front-end 
processor.  The  front-end  processor  cocwuni- 
cates  directly  with  tha  360/50  through  its 
I/O  channai,  controls  the  high-spaad  communi- 
cation lines,  and  controls  several  low-speed 
devices  such  as  a paper  tape  reader,  and 
punch  and  mark  sense  card  reader.  Other 
devices,  peripheral  to  the  360/50,^  are  con- 
trolled by  the  remote  minicomputers.  The 
total  system  permits  a wide  variety  of  I/O 
devices  which  in  turn  can  handle  data  in 
many  different  forms.  Data  can  be  input  to 
the  system  from  paper  tape,  analog  signals, 
line  drawings  (traced  on  a graphics  table) , 
pictures  (constructed  on  a graphics  terminal) , 
or  removable  disk  cartridges. 

When  the  mobile  van  is  outfitted  for 
field  experimentation,  one  of  the  Computer 
Center  NOVA's  (see  Figure  1)  is  removod  and 
installed  in  the  ven.  In  the  field,  two 
saiemometare  are  pieced  in  appropriate 
positions  and  the  ground  it  set  in  motion 
with  either  er.  explosive  charge  or  e shaker. 
Cech  seismometer  supplies  three  analog 
signals,  ons  for  esch  component  of  displace- 
ment. These  six  analog  signals  are  digitis- 


ed for  the  duration  of  the  experiment  end  the 
date  is  transforrad  to  disk.  Ths  data  can  be 
previewed  in  the  field  and  relatively  smell 
records  (*  1000  pts.  par  ohsnnsl)  can  be  ana- 
lysed se  described  in  the  previous  section. 
This  cursory  anelysie  is  vory  beneficial  for 
data  varification  and  on-aite  determination  of 
futura  seismometer  and  disturbance  placement. 


Figure  1.  Partial  Arrangement  of  UMR 
Minicomputer  Network. 

The  disk  cartridge  which  contains  the 
field  date  can  be  removed  from  tha  mobile  ver 
NOVA  end  loaded  in  a second  NOVA  located  in 
the  Computer  Center  end  networked  to  the  XIM 
360/50.  There,  the  date  is  transfarred  to  the 
360/50  disk.  At  this  point,  tha  experimenter 
osn  perform  more  sophisticated  analysis  of  the 
Oats  and  make  use  of  the  higher  performance 
peripherals.  At  present,  tha  360/50  is  used 
to  plot  the  date  and  calculation  rasulte  on  a 
Calcomp  plotter  and  punch  some  of  the  shorter 
records  on  cards  for  permanent  storage.  Long- 
er date  racords  will  be  stored  on  magnetic 
tape.  In  addition,  correlations  and  trans- 
forms can  be  conveniently  calculated  for  much 
larger  date  sets  than  can  be  processed  In  the 
minicomputer.  A future  project  will  involve 
s pattern  recognition  and  feature  extraction 
program  which  will  run  on  the  360/50  to  assist 
the  experimenter  in  his  visual  analysis  of 
data  and  calculated  results. 

Ths  total  system  permits  real-time  field 
data  acquisition,  analysis,  and  quick-look 
followed  by  central  computer  enelyeie  with 
pernenent  records  in  the  form  of  plots,  cards, 
disk  and  tape  storage. 
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SUMMARY 


INSTRUCT IOMAL  lggITg 

The  entire  data  acquisition  system  can 
ba  dismantled,  rainatallad  in  aithar  tha  cow- 
putar  cantar  or  tha  mobila  van,  and  ba  In 
operation  within  a pariod  of  fiva  hour*. 
Consequently,  thara  ia  oonsidarabla  flexibil- 
ity  in  using  tha  instrument  aa  an  acquisition 
ayatam  in  tha  field,  or  aa  an  inatructional 
and  raaaarch  inatrumant  in  tha  laboratory. 
During  tha  academic  yaar  tha  ayatam  ia  nor- 
mally uaad  for  instruction  and  processing  of 
fiald  data  acquired  during  tha  summer . 

Tha  system  is  extremely  unique  in  pro- 
viding a training  facility  for  gaophya  ;a 
students,  by  .use  of  this  inatrumant,  these  . 
students  can  aequire  a thorough  familiarity 
with  elastic  wave  methods  through  study  of 
wave  founts,  particle  motion,  correlation 
theory,  numerical  filtering,  and  deconvolu- 
tion by  manipulation  of  real  data  with  an 
imawdiata  visual  reault  of  tha  operations. 
Many  variations  of  the  filter  parameters,  for 
example,  can  be  obtained  in  a very  short 
period  of  time,  thus  giving  the  student  a 
very  sound  feeling  of  each  process,  further, 
the  student  is  only  required  to  obtain  a very 
superficial  Knowledge  of  the  computer  pro- 
grams, and  no  programming  as  ouch  is  neces- 
sary. 


The  paper  describes  tha  application  of 
minicomputers  to  tha  research  and  instruction- 
al areas  of  geophysics.  The  system  described 
here  uses  a minicomputer  in  a mobila  van  for 
tha  acquisition  and  preliminary  fiald  analysis 
of  data  produced  in  generalised  seismic  exper- 
iments. Tha  minicomputer  controls  tha  acqui- 
sition, storage  and  display  of  raw  data  and 
transformed  data.  Tha  data  is  than  manually 
transferred  to  a separate  but  similar  mini- 
computer for  loading  into  a larger  computing 
system  for  further  analysis  and  storage. 
Lastly,  the  lata  may  be  retrieved  and  display- 
ed for  classroom  use  as  a partial  substitute 
fc?  actual  field  experimentation. 
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